Avenues for Analytic exploration in Axisymmetric Spacetimes. 
Foundations and the Triad Formalism. 
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Axially symmetric spacetimes are the only models for isolated systems with continuous symmetries 
that also include dynamics. For such systems, we review the reduction of the vacuum Einstein field 
equations to their most concise form by dimensionally reducing to the three-dimensional space of 
orbits of the Killing vector, followed by a conformal rescahng. The resulting field equations can be 
written as a problem in three-dimensional gravity with a complex scalar field as source. This scalar 
field, the Ernst potential is constructed from the norm and twist of the spacelike KiUing field. In the 
case where the axial Killing vector is twist-free, we discuss the properties of the axis and simpfify 
the field equations using a triad formalism. We study two physically motivated triad choices that 
further reduce the complexity of the equations and exhibit their hierarchical structure. The first 
choice is adapted to a harmonic coordinate that asymptotes to a cylindrical radius and leads to a 
simplification of the three-dimensional Ricci tensor and the boundary conditions on the axis. We 
illustrate its properties by explicitly solving the field equations in the case of static axisymmetric 
spacetimes. The other choice of triad is based on geodesic null coordinates adapted to null infinity as 
in the Bondi formalism. We then explore the solution space of the twist-free axisymmetric vacuum 
field equations, identifying the known (unphysical) solutions together with the assumptions made in 
each case. This singles out the necessary conditions for obtaining physical solutions to the equations. 

PACS numbers: 04.20.-q, 04.20.Cv, 04.20.Jb 
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I. INTRODUCTION 

Numerical relativity has revolutionized our under- 
standing of General Relativity (GR) in the last decade, 
allowing us to study situations of high curvature and 
strongly nonlinear dynamics (see [l| for a comprehensive 
review). In particular, numerical relativity has allowed 
for the solution of the two-body problem in GR, giving a 
description of the interaction and merger of compact ob- 
jects. Successful numerical simulations of merging black 
holes have shown that these events can be well described 
by Post-Newtonian theory up until the black holes are 
quite near merger, and after merger black hole perturba- 
tion theory accurately describes the ring-down. Where 
perturbation theory fails, a simple transition between the 
regimes of the "chirp" waveform associated with Post- 
Newtonian theory and the exponential decay to a station- 
ary black hole is observed. A primary focus of current 
research is to combine these computationally expensive 
simulations with analytical approximations to create full 
inspiral-mcrgcr-ringdown gravitational waveforms 043- 
Such waveforms will serve as templates for the matched- 
filtering based signal detection methods that will be used 
in ground-based gravitational- wave detectors coming into 
operation within the next few years [ol-fl^ . 

Despite the success of perturbation and numerical 
methods in modeling binary merger waveforms, a de- 
tailed understanding of the nonlinear regime of a binary 
merger remains an open problem. It is in this stage of 
the merger that the black hole binary emits most of its 



radiated energy (see [T^ and the references therein) and 
experiences a possibly strong kick due to beamed emis- 
sion of radiation |14| - [l7| . As such, deeper analytic under- 
standing of nonlinear dynamics in GR, including better 
insights into the two-body problem, gravitational wave 
generation, and black hole formation, remains a primary 
research goal. 

The purpose of this paper is to review and expand on 
the analytic techniques involved in the study of the field 
equations in axisymmetry. Along the way we will collect 
many known and useful results, placing them into a uni- 
fied context and notation. We intend this comprehensive 
overview of the state of knowledge in the field to serve as 
a launching point for future analytic investigations and 
searches for physically relevant, exact dynamical solu- 
tions in an era where a wealth of numerical data from 
simulations is available to guide our intuition. 

We will largely restrict the discussion to the simplest 
case where there is no rotation about the axis of symme- 
try (so that the Killing vector of the symmetry is "twist- 
free"). While specializing to such a great degree docs 
limit the scope of our discussion, at least two interest- 
ing scenarios are still included in the spacetimes under 
consideration. The first is the case of a head-on merger 
of two non-rotating black holes, the simplest instance of 
the two body problem in GR. The second is the critical 
collapse of axially symmetric gravitational waves, which 
gives insights into the formation of black holes (for a re- 
view see [18[). 

Our approach to exploring the Einstein field equations 
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in this context closely follows the methods developed by 
Hoenselaers, Geroch, and Xanthopolous [l9l - l2^ . The ba- 
sic idea is to reduce the number of equations to a mini- 
mum by applying a dimensional reduction and conformal 
rescaling to the axisymmetric field equations. The result- 
ing equations are then expressed on a null basis, in the 
manner of the Ncwman-Penrosc (NP) formalism [25| but 
in only three dimensions. This triad formalism imposes 
an additional structure on the equations to be solved, 
which can lead to valuable physical insights as in the NP 
formalism. As will be illustrated in the text, the result- 
ing system of equations is simple enough to allow us to 
keep track of the assumptions made in trying to obtain a 
solution and to analyze the properties of a given solution. 
This approach may have the potential to make dynamical 
spacetime problems analytically tractable and to provide 
a consistent framework for systematically characterizing 
the results of axisymmetric numerical simulations. The 
formulation given here also has a close connection to that 
used to find solutions to the well-studied stationary ax- 
isymmetric vacuum (SAV) equations. 

To place our work in a broader context and to moti- 
vate the approach to the field equations advocated here, 
we now briefiy review the development both of the field 
of exact solutions as well as aspects of the subsequent de- 
velopment of numerical relativity. Symmetry has often 
played a primary role in arriving at a solution to the field 
equations (Ref. [26| contains a comprehensive review). 
Famous solutions such as the Schwarzschild black hole, dc 
Sitter, Anti de Sitter, and Friedmann-Robertson- Walker 
cosmological solutions all possess large numbers of sym- 
metries. Relaxing the degree of symmetries present, but 
still imposing sufficient symmetry to make headway in 
solving the field equations, leads to the study of sta- 
tionary, axisymmetric vacuum (SAV) spacetimes (cquiv- 
alcntly, spacetimes with two commuting Killing vectors) , 
which has been completely solved [ll-il U^^. It 
was shown that the task of solving the SAV equations 
can be reduced to seeking a solutions of Ernst's equa- 
tion [23, 113] on a fiat manifold. Various techniques to 
generate new solutions from known ones were developed 
in [lO, [21I, [2^ [2^ , based on examining the integral ex- 
tension (prolongation structure) of the SAV field equa- 
tions. Over the ensuing decade a variety of additional 
iques were 

ton and inverse scattering techniques [38|, and the use 
of generating functions to exponentiate the infinitesi- 
mal Hoenselaers-Kinnersley-Xanthopoulos (HKX) trans- 
formations [33! - [36j . to name a few. These techniques 
are all interrelated (39| . and each has in turn taught 
us about the structure and properties of the SAV field 
equations. They allow, for example, the generation of an 
SAV spacetime with any desired asymptotic mass- and 
current-multipole moments [4^ . Unfortunately, by their 
nature, SAV solutions cannot include gravitational radi- 
ation and tell us little about the dynamics of spacetime. 

Building on the progress made in studying in space- 



techniques were explored, including the use of harmonic 
maps [H, [3OI, Backland transformations [Sll [s^, soli- 



times with two Killing vectors, in the early and mid- 
1970 's triad methods were developed for spacetimes 
with a single symmetry, and applied to stationary 
spacetimes |4l[ and to dynamical, axisymmetric space- 
times [22i424| . At this time, however, the availability of 
increasingly powerful computers offered a promising new 
approach to obtaining solutions of the Einstein field equa- 
tions for fully generic spacetimes by numerical means. In 
the relativity community at large, the major focus of re- 
search on solving the field equations shifted from system- 
atically exploring the analytic structure to attempting 
their solution numerically. However, the numerical inte- 
gration of the field equations proved to be unexpectedly 
difficult, especially in the axisymmetric case. 

The advent of strongly hyperbolic and stable formu- 
lations of the field equations (e.g. the commonly used 
BSSN [43, and generalized harmonic 44 1 formula- 
tions) made the long term simulations of binary black 
hole simulations an exciting reality. With the steady 
progress since the breakthrough by Prctorius [491, the 
merger of compact objects has become routine [sl. l46j . 
although still computationally limited in duration and 
mass ratio. The insights afforded by these successes can 
now serve to guide research efforts aimed at obtaining an 
analytical understanding of dynamical solutions to the 
field equations. The relative simplicity of the gravita- 
tional waveforms and other observables generated during 
the highly nonlinear phase of a binary coalescence indi- 
cate that even this phase of merger could potentially be 
amenable to analytic techniques. A renewed interest in 
analytic investigations of axisymmetric spacetimes (47| 
has already led to interesting results such as the discov- 
ery and use of geometric inequalities _ lisl - fssi , studies of 
the radiation in a head-on collision 



understanding gravitational recoil [59 



models for 
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and geometri- 
es 



cal insights on gravitational radiation 

Initially, symmetries played an important role in the 
development of numerical relativity because of the great 
reduction in computational cost in axisymmetry com- 
pared to a fully 4D simulation. Some of the first suc- 
cessful work in numerical relativity was done in axisym- 
metry [6^-|68| . None of these early simulations were long- 
term stable however, due to difficulties with coordinate 
singularities at the axis of symmetry which became a 
ubiquitous issue in axisymmetric codes. Various fixes for 
the issues with the axis were implemented [69l - [7l | but 
the codes remained unstable, even when using strongly 
hyperbolic formulations of the field equations [T^I- Be- 
cause of these difficulties with axisymmetric codes, in the 
end, it was a fully 4D simulation that proved to be sta- 
ble and was able to successfully merge black holes [45j , as 
were the stable codes that followed. Successful axisym- 
metric codes capable of l ong -term evolutions have only 
recently been developed [zl, [tI] . The continued inter- 
est in axisymmetric simulations is driven mainly by the 
desire to understand the critical collapse of gravitational 
waves and in the higher accuracy and lower compu- 
tational cost of these simulations. In addition, similar 
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dimensional reductions as in the axisymmetric case are 
also used in numerical simulations of spacetimes in theo- 
ries with higher dimensions; see e.g. the mergers studied 

in [zi,!?!. 

In addition to the simplicity of the nonlinear dynamics 
observed in numerical simulations of the merger event, 
there are other tantalizing indications that the axisym- 
metric problem could be solvable analytically. The field 
equations of axisymmetric spacetimes can be written in 
terms of a generalized Ernst potential on a curved back- 
ground, given the appropriate dimensional reduction [l9} , 
discussed in detail in Sec. |lll In this formulation some of 
the techniques used to find solutions for the SAV field 
equations, such as harmonic maps, have a straightfor- 
ward generalization to the dynamic axisymmetric case. 
Viewing the numerical simulations in the context of the 
analytic techniques employed in the past may help pro- 
vide new insight into questions such as the nature of 
initial junk radiation in numerical simulations, the rea- 
sons for the robustness of certain approaches such as the 
puncture method [s^l , the nature of singularity formation 
during a collapse process, and the possible distinction be- 
tween which features of initial data contribute to the the 
mass of the final black hole and which components are ul- 
timately radiated away (similar to the way in which poles 
and scattering data can be differentiated in the nonlinear 
solution of the KdV equations [8l]|). 

The intent of this work is to provide a framework that 
could be used in future work to explore and interact with 
the results of axisymmetric simulations, drawing on the 
accumulated analytic and numerical results available for 
these spacetimes to date. We now briefiy outline the 
structure and contents of the paper. 



A. Overview of this paper 

We will begin our discussion in full generality, explic- 
itly carrying out in Sec. |IT] the series of reductions that 
ends in the field equations for vacuum, twist-free, axisym- 
metric spacetimes. Here we largely follow the discussion 
of [2^ , although in Appendix |X] we present the deriva- 
tion in the familiar notation of the 3-1-1 decomposition 
used in numerical relativity. The resulting set of equa- 
tions is equivalent to 3D GR coupled to a complex scalar 
potential £ which obeys the Ernst equation. The man- 
ifold S on which these fields are defined is obtained by 
conformally rescaling the metric on the quotient space S 
with the norm of the axial KV. The space S should be 
thought of as the physical 4D manifold M modulo the 
orbits of the Killing vector (KV) C'', or S ^ M/i''. We 
then specialize to the case of non-rotating spacetimes, 
where the field equations are equivalent to 3D GR with 
a real harmonic scalar field source that obeys the Klein- 
Gordon equation. In Sec. IIIFI we discuss general con- 
siderations regarding the existence of an axis, and note 
that the problem of divergences at the axis is in principle 
easily handled analytically. 



In Sec. lIIII we express the 3D field equations in terms of 
a triad formulation that was first developed by Hoense- 
laers f22l - [2^ , but seemingly not used by other authors (it 
should be compared to a similar triad formulation pre- 
sented by Perjes [4l[ and used in the case of stationary 
spacetimes). This formulation is derived by the projec- 
tion of the dimensionally-reduccd field equations onto a 
3D null basis, composed of two null and one orthonormal 
spatial vector. The field equations and Bianchi identities 
are then written out in full, in terms of the 3D rotation 
coefficients. Considering the success of the NP equations 
and the valuable insights they provide, this formulation 
of the axisymmetric field equations merits a more thor- 
ough investigation than it appears to have received. 

In Sec. IIVI we relate the 3D rotation coefficients and 
curvatures quantities to the familiar NP quantities on M , 
thus providing a dictionary between NP quantities and 
the quantities that arise in the triad formulation. This 
facilitates a connection to known results and an interpre- 
tation of the physical content of the triad equations. 

In the triad formulation, we have the freedom to spe- 
cialize our choice of basis vectors. In Sec. |Vl we present 
two useful choices of triad vectors and accompanying co- 
ordinates which serve to simplify the field equations. The 
first choice is, to our knowledge, new, and is analogous 
to the use of Lagrangian coordinates in fluid mechan- 
ics. In this triad choice the spatial triad leg is adapted 
to the gradient of the scalar field that encodes the dy- 
namical degree of freedom in the twist-free axisymmetric 
spacetime. By virtue of the field equations, this scalar 
field is a harmonic coordinate which asymptotically be- 
comes a cylindrical radius. This first coordinate choice 
is well-suited for analyzing the behavior of the metric 
functions and rotation coefficients on and near the axis. 
The second triad choice is inspired by the tetrad com- 
monly used in the NP formalism, where one null vector 
is taken to be geodesic and orthogonal to null hypersur- 
faces. This choice is useful in that it connects directly to 
many known solutions of the field equations, and to the 
dynamics at asymptotic null infinity, where the peeling 
property |82| - [8a | holds. 

Our purpose in Sec. IVIl is twofold. The first is to cata- 
log known axisymmetric vacuum solutions, together with 
the assumptions that lead to each solution in terms of the 
triad formalism. This isolates the conditions required for 
the spacetime to represent a physically relevant solution. 
Secondly, we provide two example derivations of (known) 
spacetimes in the context of the triad equations, to il- 
lustrate typical techniques used to find solutions in this 
formulation. While we generally do not say much about 
the extensively-studied SAV spacetimes (see e.g. [1^), in 
Sec. IVIII we discuss the equations governing SAV space- 
times in the context of our new coordinate choice from 
Sec. IV Al We conclude in Sec. IVIIII Additional useful 
results are collected in a series of appendices. 

Throughout this paper, we use geometrized units with 
G = c = \. We use Einstein summation conventions, 
with Greek indices indicating 4D coordinate indices (in 
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practice these can be taken as abstract tensor indices). 
Latin indices from the middle of the alphabet (i, j, k, . . .) 
run over 3D coordinates, and Latin indices from the be- 
ginning of the alphabet (a, 6, c, . . . ) run over 3D triad 
indices. Indices with a hat correspond to tetrad compo- 
nents of a tensor in the physical manifold Ai and run over 
1,2,3,4. Indices preceded by a comma either indicate a 
partial derivative with respect to the coordinates, as in 
f^i, or the directional derivative of a scalar quantity with 
respect to the members of a null basis, as in f a- Simi- 
larly, we use semicolons to denote covariant differentia- 
tion in the coordinate basis, as in Indices preceded 
by a bar as in Va\b indicate the intrinsic derivative on 
the triad basis. Symmetrization of indices is denoted by 
enclosing them in parenthesis, and anti-symmetrization 
by using square brackets. We use a spacctime signature 
of ( — h ++) on the 4D spacctime, and ( — h +) on the 
3D quotient space. Note that this modern convention 
differs from the signature used by many authors in the 
literature referenced here. An asterisk denotes complex 
conjugation. 



II. REDUCTION OF THE AXISYMMETRIC 
FIELD EQUATIONS 

In this section, we review a formalism for expressing 
the full four dimensional Einstein field equations in a 
simpler three dimensional form when there is a single 
continuous symmetry present in the spacctime. We then 
specialize the resulting equations to the vacuum case, and 
then to spacetimes that admit a twist-free Killing vector. 

The formalism for the dimensional reduction was pre- 
sented by Geroch for a single symmetry, and ex- 
tended by him to the case of two commuting symmetries 
[2l| in order to study SAV spacetimes. This reduction 
has been extensively used, especially in the investigation 
of stationary spacetimes [2a] . We closely follow Geroch's 
derivation and notation in what follows. We also compare 
the dimensional reduction to the familiar 3 -t- 1 decom- 
position used in numerical relativity, for which [stI . Issj 
provide excellent references. Finally, Dain's review of ax- 
isymmetric spacetimes [47| complements the discussion 
provided here and throughout this paper. 

The reduction in complexity when one studies the field 
equations for a 3D Lorentzian metric as opposed to a 4D 
metric becomes immediately apparent by counting the 
number of independent components of the Weyl tensor, 
given by N{N + 1){N + 2){N ~ 3)/12 in N dimensions. 
That is, zero independent components in 3D, and ten in 
4D. 

In the reduction to axisymmetric, vacuum spacetimes 
discussed in greater depth in Sec III Al all of the gravita- 
tional field's dynamical degrees of freedom enter as two 
scalar functions whose gradients serve as sources for the 
3D Ricci curvature. In the twist-free case, one of these 
scalars vanishes. The fact that the gravitational field is 
determined by a single remaining scalar demonstrates the 



tremendous simplification over the full 4D case with no 
symmetries present. 

The reduction proceeds in three steps. The first step is 
to derive the equations on the three manifold. Presented 
in Sec. Ill Al this process is similar to the 3-f 1 spacetime 
split familiar to numerical relativists. As a second step, 
we specialize to vacuum spacetimes. The last step of the 
reduction is a conformal rescaling, discussed in Sec. IIIDI 
which simplifies the 3D field equations further and makes 
apparent the existence of a generalized Ernst potential. 

A. The space of orbits and the general reduction of 
the field equations 

We begin by considering a 4D manifold M that ad- 
mits a metric and a Killing vector (KV) field 
Throughout this paper, we will consider to be space- 
like; however, the same formalism is easily extended to 
the case of a timelike symmetry [l^, [2^ . The KV field 
represents a continuous symmetry, and it defines a set 
of integral curves called the orbits of Motion along 
these orbits leaves the spacetime invariant and preserves 
the metric. This means that tensor fields on Ai have 
vanishing Lie derivative along For the case of the 
metric tensor /^^.g^jy = leads to the Killing equation, 

V(^e.) = o. (2.1) 

Intuitively, we see that one of the dimensions of Ai is 
redundant, and so we would like to reduce the study of 
this spacetime to the study of some 3D space. 

Naively, one would think of considering dynamics in 
Ai only on surfaces to which is orthogonal. In prac- 
tice however, is only orthogonal to such a foliation of 
submanifolds of Ai if its twist oj'', given by 

- e^^.p^r VT , (2.2) 

vanishes. When = 0, the KV S,^ points in the same 
direction as the gradient of some scalar function (/) on Ai, 
but this is not true in general jS^] . Instead of considering 
some hypersurface in Ai , we consider a new space, which 
we call S following Geroch . The space S is defined as 
the collection of orbits of in A^; it is a 3D space that 
can be shown to posses all the properties of a manifold. 
The space S can be represented as a surface in Ai only 
if = 0. Figure [T] provides an illustration of the case of 
a twist-free symmetry with closed orbits. 

We denote with an over-bar tensor fields on S. These 
fields are orthogonal to the KV on all their indices, e.g. 
T^^^ = T^^„ = 0. A metric /i^^ on S can be defined 
by "subtracting" the exterior product of two unit vectors 
pointing in the direction of the KV from the metric g^,^. 
The resultant metric on S is 

hf„y = - A^^^^^i. . (2.3) 

Note that £,^h^i, ~ and the Lie derivative of h^i, along 
vanishes. The function A that appears in Eq. (|2.3p is 
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also a subspace of M. The fact that the orbits are closed 
implies that a set of fixed points, namely the axis, must exist 
if the spacetime is asymptotically flat. 

the norm of the spacelike KV, 

= A > , (2.4) 

and will play a key role in the reduction that follows. 

By raising an index on /i^^ using g^'^ , we can define 
a projection operator h" which projects 4D fields onto 
S. Arbitrary tensor fields can be projected into S by 
contracting all of their indices onto the projector, 

= and f^p^K'h-'^T^,, (2.5) 

and similarly for tensors of arbitrary rank. We also define 
the operator Da by contracting the usual 4D covariant 
derivative of a tensor field with the projector on all its 
indices, 

Dafp^ = h^h'phP{VM. (2.6) 

It can be shown that the operator Da obeys all the usual 
axioms associated with the unique covariant derivative 
operator on a manifold with metric /i^^ [20j . 

Given the metric h^^, on S and a compatible covari- 
ant derivative, we can compute the Riemann tensor on 
S, and relate it to the 4D Riemann tensor and the KV 



^f^. In doing so, the 4D field equations will be expressed 
entirely in terms of quantities on S. This projection of 
the 4D field equations is achieved by writing out Gauss- 
Codazzi equations generalized to the case of a timelike 
quotient space. This calculation, although computation- 
ally intensive, is only a slight modification of the standard 
techniques of the 3-1-1 split often used in numerical rela- 
tivity and is detailed in AppendixlX] Here, we summarize 
the key results that will be used later in the text. 

The contracted Gauss equation expresses the 3D Ricci 
curvature Rap on S in terms of the Ricci tensor Rf^i, on 
the manifold A4 , derivatives of the norm A of the K V and 
its twist oj'' as 

Rap ^K'h'pR,.. + ^DaDpX - ^Da\Dp\ 

- ^ [hapLO^Lo'' - WaW/j) • (2.7) 

Since 5 is a 3D manifold all the curvature information 
on S is contained in the Ricci tensor Rap associated with 
h^ui with the remaining geometric content of given by 
the magnitude A and twist w'^ of Note that Eq. ((2J)) 
has the same form as the Einstein field equations on the 
three manifold S with additional source terms on the 
right hand side; in the case where there are matter fields, 
we would re-express R^^ in terms of the stress energy 
tensor T^i,. We are primarily interested in the vacuum 
field equations, in which case i?^^ = and the geometry 
on the three manifold is entirely sourced by A and u)^. 
As such, we need equations governing the evolution of A 
and uj^ in order to complete our reduction of the field 
equations. 

This second set of equations is analogous to the Co- 
dazzi equations [87| . since they are derived by applying 
the Ricci identity to the unit vector tangent to the KV. 
They are detailed in Appendix lA 21 The resulting equa- 
tion governing A is 

^'A ^^DaXD^X - jio^uj^ - 2i?^,^^r , (2.8) 

where the 3D wave operator is defined using D^ = 
D^Da- The twist obeys the equations 

D"LOa=-^UaD^X, (2.9) 
D[ai^^]=-eaPpa^''R"rC- (2.10) 

Together, Eqs. (|2.7p -- (|2.10p can be solved on S for hap, A 
and ujfj,. We can then find an expression for the KV S,^ 
using the identity, derived in Appendix I A 21 

v^e. = ^e^.p^e^c^" - ^e[^v,]A (2.11) 

together with the fact that ^'^ft.^i/ = 0. With the KV and 
/ipi/, we can finally reconstruct the full 4D metric g^^i, on 
Ai, completing the solution of the field equations. 

The field equations on S are greatly simplified com- 
pared to the full Einstein field equations, but they are 
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still formidable. As such, we will make a series of further 
specializations with the aim of rendering them tractable. 
In the past, the assumption of a second, timelike symme- 
try has resulted in the SAV equations and their solution. 
We will briefly discuss the SAV equations in Sec. I VIII 
in the context of a convenient coordinate system we in- 
troduce in Sec. IV Al Since our purpose is to pursue new 
solutions, outside of Sec. I VII] we will not assume any fur- 
ther symmetries. Instead, we give the reductions of the 
field equations in the case of vacuum, and then twist-free 
spacetimes in the sections that follow. 



B. Coordinates adapted to the symmetry 



This decomposition of the 4D metric and its inverse in 
terms of the 3D metric and the KV should be compared 
to the analogous decompositions of the 4D metric into 
a spatial metric, lapse, and shift vector in a 3 -I- 1 split. 



e.g. as found in [87|. For the remainder of this text. 



we will use coordinates adapted to the Killing symme- 
try, so that the decompositions (j2.13p and (|2.14p hold. 
The most useful consequence of this choice is that all of 
the information contained in quantities projected onto S 
is contained in the components on the coordinate basis 
. As such, we will write projected four dimensional in- 
dices a, /3, . . . as Latin three dimensional indices such as 
I, J, fc, . . . which run over coordinates on S. 



In this section we detail the consequences of using 
a coordinate system adapted to the Killing symmetry. 
For a spacetime admitting a KV there exists coordi- 
nates x^^ = (x* on such that ~ S'^, where 
is a coordinate that does not appear in the metric, 
C^g^^ = dg^^/d(j) = [M [93- To find the form of the 
metric g^^i, in coordinates adapted to an axial KV, we 
first note that 



(2.12) 



which also implies that = A. We denote the remaining 
covariant components of by -Bi, so that = (i?i, A). 
Since fully projected quantities on S arc orthogonal to 
e.g. VaSf^ = = 0, the components of projected 
tensors vanish, and the remaining components of h^^ are 
the 3x3 block of components htj . Using this in Eq. 
the metric 17^^ takes a simple form 



h,j + X- 
B, 



'B,B., 



B, 
X 



(2.13) 



C. The vacuum field equations 

We now consider the case of vacuum 4D spacetimes. 
This sets the 4D Ricci tensor to zero in the equations 
derived in Sec. Ill Al Importantly, we see from Eq. (|2.10p 
that the curl of the twist vector vanishes. We can thus 
define a twist potential w such that 



(2.17) 



From Eqs. (f^ . (f^ . and (PT7)) . rccalhng that we may 
use 3D indices for quantities projected onto S, we have 
as our field equations 



1 



1 



D^X = —D'XD^X - -D'ujD.uj , 
2A A 

D^uj ^ —D,ujD'X, 
2A 

+ ^D,D,X-^D,XD,X. 



(2.18) 



Denoting the inverse of hij by h^^ and using it to raise 
and lower 3D indices, we can define i?* = h^-' Bj and 



B' 



h^^ BiBj. This allows us to write the inverse of the 



metric (|2.13p as 



-X-^Bj (A-S2)-i 



The determinant of can be expressed as 



det = g = X det hi 



Xh. 



(2.14) 



(2.15) 



Finally, in this basis the relationship between twist of 
the KV and Bi can be found by defining the projected 
antisymmetric tensor ea/3-y = tafjjfj.S.'^ /VX- Using the 
definition of the twist (|2.2p and projecting onto S we 
have 



LUi = y/Xe^^'^DjBk- 



(2.16) 



from which we can see that if Bi vanishes, so does the 
twist. 



D. The conformally rescaled equations and the 
Ernst potential 

A further simplification to the reduced field equa- 
tions (|2.18p can be obtained by conformally rescaling the 
metric hij. We define ft-^^ to be 



(2.19) 



and investigate the conformally rescaled 3D manifold 
which we will call S. The vacuum field equations (j2.18p 
can now be rewritten in terms of hij , bearing in mind that 
the Christoffel symbols associated with the two metrics 
are related by 

ni = + ^ {S]X^k + SIX^, - h,kh''X,i) . (2.20) 

The wave operator, associated with h^p is related to 
Z)2 as 



(2.21) 
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and further, D'fDJ ^ XD'fDif. Substituting these 
identities into Eqs. (|2.18p . the field equations can be ex- 
pressed using the metric hij [13, 

A A 



D^Lo = -ujiD'X, 
X 

Rf = ^ [D.ujDjLj + D.XDjX] 



(2.22) 



The symbol R^J' denotes the Ricci curvature of the 
rescaled three manifold with metric hij. There is ad- 
ditional structure in these equations which can be made 
more apparent ^by introducing the complex Ernst poten- 
tial £ ^ X + iuj [2i\. In terms of this potential, Eqs. ((Z!22)) 
become 



D^E = 



2D,£D'£ 
{£ + £*) ' 



Rf? = 



2D(i£Dj)£* 
{£ + £*Y 



(2.23) 



It is important to note that the Ernst potential usually 
discussed in the context of stationary spacetimes is based 
on the norm and twist of a timelike KV, rather than the 
spacelike KV as discussed in this section. This results 
in some sign differences in various definitions, c.f. the 
relevant chapters of [2^. The relationship between the 
Ernst potential defined here and the Ernst potential used 
in conjunction with SAV spacetimes is explained further 
in Sec. IVlIl 



E. 



Reduction to the case of twist-free Killing 
vectors 



The axisymmetric field equations (j2.23p , though much 
simplified from their full 4D form, remain intractable. 
For the remainder of this paper, we restrict our explo- 
ration to the situation depicted in Fig. [TJ where is 
hypersurface orthogonal, so that w = 0. In doing so we 
eliminate the possibility of the study of rotating axisym- 
metric spacetimes, but we benefit from further simplifi- 
cations to the field equations. A number of physically 
interesting dynamical spacetime solutions are twist free, 
these including the head-on collision of black holes and 
non-spinning, axisymmetric critical collapse. 

The twist-free assumption reduces the problem of find- 
ing solutions to the field equations to the study of a har- 
monic scalar ip on the three manifold S, where we define 
■0 via 



A = e"^ 



The field equations (j2.23p become 



= 0, 



(2.24) 



(2.25) 



and the Ricci scalar associated with the three met- 
ric, which we denote as i?, is given by the contraction 

of ^IM. 



where we used semicolons in place of Di to condense the 
notation for the covariant derivatives. The scalar R is 
the only nonzero eigenvalue of and corresponds to 
the eigenvector ip.^i. 

Some general properties of gravity in 3D are discussed 
in [9l|. In particular, since the 3D gravitational field has 
no dynamics, due to the vanishing of the Weyl tensor, the 
only dynamical degree of freedom in the problem is the 
scalar tp. This reduced number of variables drastically 
simplifies the calculations. In the sections that follow we 
present a systematic way of analyzing Eqs. (j2.25p using 
a triad formalism, without immediately specializing to 
any given coordinate system. The fact that ^ is har- 
monic makes it a convenient choice of coordinate on S 
that in addition greatly simplifies the components of the 
Ricci tensor i?fj^. The full implication of choosing -0 as 
a coordinate, as well as another gauge choice adapted 
to geodesic null coordinates on S, are discussed in the 
sections that follow. 



F. The axis 

All of the previous results in this section hold for KVs 
with generic orbits. Here, we review some additional re- 
sults which apply if the orbits are closed, as in the case 
of axisymmctry. Motion along the orbits of a KV maps 
the spacetime onto itself, and by the definition of the KV 
this map preserves the metric. This map may have fixed 
points, where it is simply the identity operator, and these 
fixed points comprise the axis of the spacetime. Much is 
known about the axis of an axisymmetric spacetime, see 
e.g. [92l - i94j . A key result due to Carter is that any 
vacuum spacetime with a KV that has closed orbits and 
is asymptotically flat admits fixed points and therefore 
isolated systems which possess an axial KV will have 
an axis. 

This axis is 2D and timelike [o^l, and will be denoted 
W2- On the axis the magnitude of the axial KV vanishes. 



W2 



0. 



(2.27) 



R = 20;,^/>'' , 



(2.26) 



Note that the derivative of the KV, cannot vanish on 
the axis, or else would vanish everywhere (see e.g. [s^ 
for further discussion). 

When the axis is free of singularities, a condition 
known as elementary flatness holds in a neighborhood 
of the axis. This condition expresses the fact that in 
the local Lorentz frame of a small neighborhood about a 
point on W2, we can make a loop around the axis, and 
the circumference of this loop must be equal to 2-k times 
its radius. If this is not true, then there is a conical sin- 
gularity in this small neighborhood, and traversing the 
circle around them results in a deflcit (or surplus) an- 
gle. One way to express elementary flatness is to find a 
set of coordinates in which the line element has the form 
= gppdp^ + XdcfP' near the axis, holding the third spa- 
tial coordinate fixed. Dividing the proper length around 
a circle by 2tt times the proper distance to the axis yields 
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a constant Kj^ which, if different from unity, gives a mea- 
sure of the deficit angle [95| . 



Kr 



lim 



27r /p'' ^5;;^ dp 



(2.28) 



A coordinate invariant form of this same condition that 
is more useful from our perspective was given by Mars 
and Senovilla I93i. 



hm^ 

A-^o 4A 



1. 



(2.29) 



We derive this result using a specific coordinate system 
in Sec. IV Al Expressing Eq. p.29p in terms of ip and the 



conformal three metric ha we have 



lim e 

A->0 



4ip 



(2.30) 



Equation (j2.30p provides explicit boundary conditions for 
quantities on S as the axis A = is approached, if we wish 
our axis to be free of conical singularities. 



III. THE TWIST-FREE FIELD EQUATIONS 
EXPRESSED USING A TRIAD FORMALISM 

To explore the field equations on the 3D manifold S, 
we employ a triad formalism in which we choose a basis 
for the tangent bundle before further selecting coordi- 
nates on the manifold. In this section we follow Hocn- 
saelers [2^ - |2^ in writing out the 3D field equations (|2.25p 
and Bianchi identities on 5 in a manner similar to the 
Newman- Penrose (NP) equations [l^lsBl- This form of 
the equations is particularly convenient for the study of 
the exact solutions of the field equations, since it makes 
manifest what the various possible assumptions and sim- 
plifications might be for special and physically interesting 
cases. The procedure is to define a null (or orthonor- 
mal) triad and write out in full the field equations ex- 
pressed in this basis (we note that a similar formalism 
was developed by Perjes in [4l| in the context of station- 
ary spacetimes, using a complex triad). Our approach 
largely follows the conventions for the tetrad formalism 
used in Chandrasekhar's text [9^, which also gives gen- 
eral background on the technique. 

We begin by selecting a triad basis 

C = (r,n\cO, (3.1) 
so that the metric expressed on this triad basis 

Vab = h,,ClCi (3.2) 

contains only constant coefficients. A null triad choice 
that is particularly useful is one for which 

lil'' = riin^ = lid = riic' = and qc* = —Un^ ~ 1 . 

(3.3) 



The orientation of the triad is fixed by the equations 



f-ijk 



(3.4) 



Given the normalization in Eq. p.3|) . the metric on the 
coordinate basis is expressed in terms of the triad vectors 
as 



liTlj 



Uilj 



C{Cj 



(3.5) 



The fundamental variables in a triad formalism arc the 
Ricci rotation coefficients jabc, which record how the ba- 
sis vectors change as we traverse the manifold. They are 
defined by 



Jabc 



Caj;kCl,Cc 



(3.6) 



The rotation coefficients are antisymmetric in the first 
two indices jabc = J[ab]c (note our ordering of indices 
induces a sign change from Chandrasekhar's definition 
[96^). In 3D there are nine independent real rotation co- 
efficients, as opposed to the 24 real rotation coefficients 
that exist in 4D. We adopt the following naming conven- 
tion first introduced in 1221 , 



a 


= 7121 


— ll-jJl , 




= 7311 




7 


= 7231 


— Tli-jC l'^ , 


5 


= 7122 




e 


= 7312 


— Cj-j/ Hp , 


c 


= 7232 




V 


= 7123 


— ^v.j'^ ^ 1 


e 


= 7313 




L 


= 7233 


Ti^JjC CP • 









(3.7) 

The projection of the 3D Ricci tensor Rij (we drop the 
superscript 3D from here on) onto this basis gives us six 
curvature scalars, which we denote 



i?ll 




4>i 


~ ^12 






-R13 


' " ^ ^ 1 


02 


= ^22 


= Rijn^n-' 




-^23 




00 


= ^33 


— t^ij c c 


(3 



The Ricci scalar is given by^ 

R = R\ = ( 



204 



(3.9) 



The field equations describe how the rotation coeffi- 
cients listed in Eqs. (j3.7[) change in a particular basis 
direction to ensure that Eqs. (|2.25p are satisfied. The 
change along a basis direction is a directional derivative 
given by 



^ai---an,b — (^i---Zti Cai ^ ' ' ' Can ^ )',jCb^ 



(3.10) 



^ Note that here we are using a different definition of </>o and 04 
than [23| . If we denote the scalars of |22|| with a superscript 
H, the relationship between the two conventions is such that 
p4 and 2(p^ = 0o- 
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The basis-dependent directional derivative can be related 
to the intrinsic covariant derivative of a tensor projected 
onto the triad basis, 

Va,-a^\b = {V^,...^J,,Ca:' ■■■CaJ-Cb' , (3.11) 

by taking into account the manner in which the basis 
itself changes. Using Eqs. (|3^ . (|3A0|) . and (|3?TT1) one 
can show that the relationship between the directional 
and intrinsic derivatives is 

Ki---a„|h = Vai...a„,b + aib^c- ■ -a^ H h ar^b^a^- ■ -c ■ 

(3.12) 

Recall that triad indices are raised using the constant 
metric 77°'' defined by r]°'^riij(. = (5°, which has the same 



component form as rjab- It is important to note that in 
a triad formalism, the directional derivatives of a scalar 
function do not commute, while intrinsic derivatives do. 
The commutation relations for directional derivatives are 

f.[ab] = f.ml[a"'b] ■ (3.13) 

Using the relationship between intrinsic and directional 
derivatives, we now express the field equations on the 
triad basis in terms of the rotation coefficients The 
Ricci tensor obeys 

-^ab — la m,b I ab,m ] nlniab la ibnm- 

(3.14) 

Writing out the field equations p.l4p in full leads to 



(e-7),3H 




-02 


= m - 2/3C + 7^ + 4- 


- e{d + 2l) + 00 , 










- L{e - ?/) + 01 , 




C,3 


- ^2 


= t(i - ,5) + C(e + 2?7 - 


7) + 02 , 




a, 3 


- '7,1 


= -aij + 7]) + p{S - l] 


) + 61(7 - 77) + 03 , 


{a + 9)^2 - 




- <=,3 


= -7(r/ - e) - d{9 + 2a) - e{e + rj) -dL + (j>4 




^,1 




= 9{e-a)+l3{j + 2r] 


- e) + </'5 , 


M +^'.2 - 


(e^ 


-7),3 


= ai + 7^ - (56* - , 




(e + '?),! - 




- 


= a{-/ + rj) + l3{6 + l) - 




Ci + - 


(7^ 


-v):2 


= C{a + 0) + S{e + ?7) ^ 


- '-(7 + v) ■ 



(3.15) 
(3.16) 
(3.17) 
(3.18) 
(3.19) 
(3.20) 
(3.21) 
(3.22) 
(3.23) 

In the twist-free case, the curvature scalars 0j appearing in the above expressions are obtained from the Ricci tensor 
Rab computed by projecting Eqs. (|2.25|) onto the triad, 

Rab = 2^^a^,b. (3.24) 

Of the above set of nine field equations there are six equations that contain Ricci curvature components and three 
that do not. These three equations constitute the 3D version of the eliminant relations (see Chandrasekhar [1^). As 
expected, there are fewer equations on S than in the 4D case (9 here versus 36 equations in 4D). 
The three, 3D Bianchi identities 

Ra\b - \R,a + IbamR^"^ + Ib^ ^Ra = (3.25) 



are written as^ 



'o),i + (05),2 - (<^3),3 - ^(-^0 + 04) - (t + 2(5)05 + (?/ + 7 - 2e)03 " /301 , (3.26) 
'o),2 + (-^2),! - (0l),3 = -'•(00 + 04) + (2a + 0)02 -f (27 - 6 - ?/)0i + C03 , (3.27) 
^3),2 = -0i(20 + a) + /?02 + (e " 7)(</'4 + M + {& + 2003 - C05 . (3.28) 
I 

In terms of the rotation coefficients, the commutation f ^i — f j^ = (7 + v)!,! — + Of .3, 

relations are f _ j = _^f^^ + (, + + . (3 29) 

/.21 — /.12 = <5/.i + a/_2 + (7 + e)/.3, and must be used whenever interchanging the order of 
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directional derivatives. Finally, it is useful to note that 
the operator D^f = f J'^ can be expressed as 



-2/.12 + /,33 + 2S)f,i + df,2 - 2e/,3 

= -2/.21 + /.33 - ^li + (2a + e)f,2 + 27/,3 . 

(3.30) 

This concludes the general triad formulation using the 
rotation coefficients as fundamental variables. The equa- 
tions given are valid for both twisting and twist-free 
spacetimes, with the only difference being the complexity 
of the 3D Ricci tensor. These equations can be simplified 
to a great degree by a judicious choice of triad. We ex- 
plore two especially useful triad choices in Sec.lVH where 
we also specialize to the case of twist-free spacetimes. 



IV. RELATING PHYSICAL 4D QUANTITIES 
TO COMPUTED 3D QUANTITIES 

In this section we provide the explicit correspondence 
between NP quantities on the physical 4D spacetime Ai 
and the computationally concise quantities on the confor- 
mal manifold S. Knowledge of this correspondence is use- 
ful for various reasons: (i) Initial conditions for integrat- 
ing the much simpler 3D field equations (|3.15p - (|3.23p are 
most readily specified on M. (ii) The boundary condi- 
tions on the axis, discussed in Sec. Ill Fl and Appendix iDl 
require information about the smoothness of the physi- 
cal quantities on Ai, since the 3D conformal metric hij 
is singular on the axis, (iii) Searching for solutions to 
the field equations involves making choices of the triad 
and the gauge, and having a direct translation of the as- 
sumptions made in 3D to the implications for the physi- 
cal quantities is advantageous. This relationship between 
specializations in 3D and 4D also identifies the conditions 
on the 3D quantities corresponding to known solutions. 

To exhibit the correspondence, we first note that in 
the twist free case, the metric decomposition of Sec. Ill Bl 
simplifies to the case where Bi = and the 4D metric 
5^1/ can thus be expressed as 











(4.1) 



The metric p^^, which is conformal to the physical met- 
ric, provides a useful intermediate step for the calcula- 
tions that follow. 



^ Note that Eqs. 1 13. 261 1-1 13. 281 1 differ from those derived by Hoense- 
laers in l22l . Eq. {3.5c), with respect to the sign in front of his 
d operator. In addition to this difference and the difference in 
notation the equations presented in this section differ from those 
of Hoenselaers in that our triad has a different normalization, 
and so some of the signs differ. Specifically, every factor of I' 
and derivative in the P direction receives a sign change, which 
changes the signs in front of many of the rotation coefficients 
and some of the curvature scalars. With these considerations, 
the two sets of equations are identical. 



A. Spin coefficients 

We now define the relationship between the NP spin 
coefficients on Ai and the rotation coefficients defined in 
Eqs. (|3.7|) . There is some freedom in the choice of tetrad 
as we go between the 4D and 3D manifolds, which we fix 
by choosing the tetrad so that the directions of all the 
null basis vectors coincide, and so that the parametriza- 
tion of the out-going null vectors are the same. In order 
to avoid confusion, all quantities on A4 such as spin co- 
efficients, K, e and Wcyl scalars, are given with a hat 
(?). Quantities associated with the conformally rescaled 
4-mctric g^j^ arc all indicated with a tilde (T), and 3D 
quantities will remain unadorned. 

The standard complex null tetrad on Ai is 



(4.2) 



with the non-zero metric components being t};,. 



Vnl = 

-1. Now consider another tetrad 



constructed by augmenting the triad (j3.ip with the vector 

= e~'^^5'^ which has the same direction as the KV 
to yield the tetrad 



c = (^^n^c^d^), 



(4.3) 



where we have omitted the tilde's to emphasize that this 
tetrad is built from the same triad vectors that we use 
on S (although strictly speaking they are the lift of these 
vectors onto a conformal 4D space). It can be verified 
directly using Eqs. (|3.5p and (|4.ip that the conformal 
metric g^i, can be expressed as 



(4.4) 



where the covector is c?^ = e^^Sf^. 

To find the relationship between the NP spin coef- 
ficients on Ai and the rotation coefficients on S, we 
first calculate the rotation coefficients associated with the 
conformally related metric g^i, on the basis in Eq. (j4.3p . 
To do this expediently we introduce the quantities Xatc 
that are defined as [90| 



76a 



7bc 



(^7.^ - Cbf}n)C^C2. (4.5) 



and are antisymmetric in the first and third indices. The 
major advantage of working with the quantities Xatc is 
that they can be computed using coordinate derivatives 
rather than covariant derivatives. This property makes 
the comparison between quantities defined on different 
metrics given the same coordinate choice easy. Given a 
set of Xabc's the rotation coefficients can be constructed 
using the relation 



"fabc 



'{Xabc + Xcab ~ Xbca) 



(4.6) 



The 24 rotation coefficients associated with the con- 
formal metric g^^ can be related to the nine rotation co- 
efficients associated with hab by noting that Xabc = Xabc 
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when a,b,c run over 1,2,3. The remaining 15 rotation 
coefScients can be subdivided into nine coefficients of the 
form 7a4b, three ^abi coefficients, and three 7^44 coeffi- 
cients . From the definitions in Eq. (|4.5p and the vector 
rf'^ , it is straightforward to verify that AaM = Xaib = 0, 
and so the 12 coefficients jaib and ^abi vanish. There are 
then only three non-zero rotation coefficients, 



7a44 



-A 



a44 



(4.7) 



in addition to those in Eq. p.7p . 

Given the rotation coefficients associated with the aug- 
mented tetrad in Eq. (|4.3p . the spin coefficients associ- 
ated with the physical space tetrad in Eq. (|4.2p can be 
obtained from a transformation of the form 



Cap = p!mb^^■ 



(4.8) 



Specifically, = e ^'^Q^ and the nonzero components 
of Qa 



V2' 



e3V 



V2 



(VA3 + e) : 



(2^,3 + V), P- 



V2 



(7 -"0,3), 



2^/2 



A=-e2'/'(,-2^,2), /i=2^e2^ 



--(Je^"^. (4.12) 



The identifications in Eqs. (|4.12p gives us the benefit of all 
the usual intuition regarding the spin coefficients in the 
4D spacetime when computing quantities on the manifold 
S. We will explore these relationships and their physical 
implications more fully in Sec. I VII when we review the 
exact solutions to the field equations. 



2V' 



Q\ = -Q 



(4.9) 



Note that the fact that Q\ = 1 ensures that the 
parametrization of outgoing null vector V' on the three 
manifold coincides with the associated vector on the 4D 
spacetime. The vectors on M are then given in terms of 
the tetrad (gSl) by 



(4.10) 



By repeatedly using the definition (|4.5p on the different 
tetrads, we find that the Xabc functions associated with 
the physical tetrad (|4.2p [and thus the rotation coeffi- 
cients via Eq. (|4.6p ] are related to those on the augmented 
tetrad given in Eq. (|4.3p by 



QlQl ma{Ph,c-flbciP').a 



(4.11) 



where the constant metric fjab has the non-zero compo- 
nents, ^12 = ry2i = -1 and 7733 = fj^^ = I. 

Since the Pg°'s are functions only of ip all the physical 
rotation coefficients reconstructed using Eq. (j4.6p given 
Eq. (|4.11[) can be written in terms of the nine rotation co- 
efficients on the triad basis, the three directional deriva- 
tives of the scalar function ip, and functions of ip itself. It 
can also be observed that all the physical rotation coeffi- 
cients expressed on the basis in Eq. (|4.10p are real. The 
physical spin coefficients using the NP naming convention 
[96| , when expressed in terms of the rotation coefficients 
defined on S are 



1 



2' 



(2^,1+^): 



B. Curvature and Weyl scalars 

The second set of quantities that are useful for explor- 
ing the physical content of spacetime, such as gravita- 
tional radiation, are the Weyl scalars. In this section 
we will show that they have a particularly simple rep- 
resentation in terms of the 3D rotation coefficients and 
directional derivatives of tp. 

The fact that the Weyl tensor is conformally invariant 
implies that on the coordinate basis C"/3js = C°' p-yS- 
Lowering the index a, expressing the tensor on the tetrad 
basis in Eq. (j4.10p . and subsequently using Eq. (j4.8p to 
express it on the augmented basis in Eq. (|4.3p , we obtain 
an expression for the physical Weyl tensor in terms of the 
Weyl tensor on the augmented basis. 



CabcA = ^-^''CabMQlQlQ'\. 



(4.13) 



The quantity Cabcd is readily computed in terms of the 
rotation coefficients and directional derivatives of "0 on S 
from the standard expression for the Riemann tensor [9^ , 
which in vacuum is identical to the Weyl tensor: 



Rabcd = labcA — labd.c + f]^ ^ lbaf{lcgd — Idgc) 



fj^^ilfaclbgd - Ifadlbgc) ■ 



(4.14) 



Writing out Eqs. (|4.14[) in full, making use of the defini- 
tions of (f>i given in Eqs. (|3.8|) and (|3.24p . and substituting 
in the field equations p.l5p - p.23p wherever necessary 
yields the following expressions for the Weyl scalars on 
the physical manifold: 

^"0 = ^1313 = {ai^,i + 3(0,1)^ + i^M + P^p,3) , 
^-i = C1213 = ^ (^4 (SV'.a - 7) + i\3i + , 



c 



1342 



(7/>,l (2^,2-0+^V^.2 + 2(V',3)'+^. 



33 j 



12 



"^4 



Ci2 



V2 



- ■0,2 (3-0,3 + e) - 0,32) , 



C2424 = -e"^ - 3(^,2)' - 0,22 + CV',3) 



(4.15) 



It is important to note that the assumption of twist-free 
axisymmetry greatly decreases the number of indepen- 
dent functions to be considered: the NP spin coefficients 
and Weyl scalars which in general arc complex are all real 
in the twist-free case, effectively cutting the problem of 
finding solutions in half. Further simplifications can be 
achieved with specific gauge and tetrad choices. 



03 = 02 = 01 = and 00 = i?. This greatly simplifies 
the Bianchi identities, which are 



R 



= 29, 



R 



= -2t, 



R 



= 2(e-7), (5.2) 



and which gives the rotation coefficients appearing in 
Eq. (|5.2p the interpretation of being proportional to the 
rate of change of In R in a particular direction. 

Because i? is a scalar, the curl of its gradient, 
€°'^'^R\}jc = 0, must vanish. Equivalently, the commutator 
equations p.29p with f = R must hold. This augments 
the field equations with the following three equations. 



V. TWO TRIAD CHOICES 

In this section we discuss the implications of two 
physically-motivated triad choices which further sim- 
pfify Eqs. (|3?T5|) " (|3?23)) and the Bianchi identities dSSS])- 
(|3.28p . The first choice is to use ?/> as a coordinate and to 
associate the triad direction Ca with its gradient. This 
choice greatly simplifies the Ricci tensor on the three 
manifold and is suited to applying the boundary con- 
dition on the axis. The second is to use geodesic null co- 
ordinates. This allows us to make direct contact with the 
Bondi formalism and thus the emitted radiation reaching 
future null infinity 1+ in asymptotically flat spacctimcs. 



A. Choosing ?/j as a coordinate 



The field equations p.25|) describe a gravitational field 
on a three manifold sourced by a harmonic scalar field ip 
which obeys D^ip = 0. In 4D gravity, harmonic coordi- 
nates have been successfully employed, e.g. for proving 
the wcU-posedness of the Cauchy problem for the Ein- 
stein equations [s^ [o^l ■ The usefulness of harmonic co- 
ordinates in 4D, together with the fact that the 3D Ricci 
tensor greatly simplifies if ip is chosen as a coordinate 
leads us to investigate this gauge choice further. 

We now specialize our triad so that points in the 
same direction as the gradient of ip- The normalization 
condition Cqc" = 1 implies that 



(5.1) 



where R = 20,^0''° is the 3D Ricci scalar defined 
in ()2.26|) . Note that the sign of R determines whether -0 q 
is timelike, spacelike, or null. For Schwarzschild, R > 0, 
and so we might expect this to be true of a physically 
reasonable spacetime, especially one that settles down 
to Schwarzschild after some dynamical evolution, and as 
such we will assume that tjj^a is spacelike. 

Given the definition of Ca in Eq. (|5.ip we can express 
the Ricci tensor (|2.25|) as R 
curvature scalars defined 




in Eqs. 



RciCj, and so the six 



(e -7), 1-^,3 - 
(e - 7), 2 + '.,3 - 



t(7- 



7=0, 
- V) = 0, 



(5.3) 



The fact that Ca points along the gradient of a scalar 
places additional conditions on the rotation coefficients. 
To see this, we compute the intrinsic derivative of Ca and 
express the result on the triad basis to obtain 



Ca\b 



-7 -C 





(5.4) 



Now, noting that c^ib + ^Ca{lnR),i, = y/2/Rip\abi and 
using the directional derivatives of R computed in 
we have 



[3 e 

-1 



(5.5) 



— t e — 7 



However since ip^a is a gradient, this matrix should be 
symmetric. Thus 7 = — e. Further, wc note that D'^%j: = 
is automatically satisfied. 

The Bianchi identities (|5.2|) . in addition to the field 
equations p.l5p - p.23p . allow us to find a particularly 
simple expression for the wave operator of In R, 



D^{\nR) = 2{R^2e^ ~2CP) . 



(5.6) 



It is interesting to note that if ip is chosen as a co- 
ordinate and the tetrad leg Ca is fixed using (|5.1|) . then 
the directional derivatives of -0 that enter into the 4D ex- 
pressions for the NP scalars become particularly simple. 
Explicitly tp^i = t/)_2 = and -0,3 = \fRj2. This implies 
that the expressions for the Weyl scalars (|4.15p become 




*4 



2^' 



(5.7) 
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The rotation cocfflciciits that enter these expressions are 
the same rotation coefBcients that appear in the second 
derivative of -0 expressed on the triad basis, Eq. (|5.5|) . 
This underscores the fact that the scalar tp sources the 
gravitational field. Another important consequence of 
Eqs. (|5.7p is that for this tetrad choice, if l"^ is geodesic, 
i.e. /3 = 0, then the geodesic is a principle null geodesic 
of the spacetime, = 0. 

Thus far the other triad vectors are unspecified, except 
that they are null and orthogonal to Cq. With Ca fixed, we 
still have freedom to boost along 1°-. The equivalent of the 
Lorentz transformations for the 3D triad are discussed 
fully in Appendix |B] Here we consider the effect of a 
boost of the form 



(5.8) 



Using the definitions in Eqs. (j3.7|) . we find that under 
such a boost, six of the coefficients are simply multiplied 
by factors of A, while three have nontrivial transforms. 



a = Aa — A,i fj 



At] - A^3 
A 



~ AS - A^2 



The full transforms are given in Eq. (jB2[) : interestingly, 
the above coefficients with a nontrivial transform do 
not enter into the expressions for the Weyl scalars in 
Eqs. (l5Jll . 

We can always use our boost freedom to set at least 
one of a, 77, or ^ to zero. Note that if a boost exists that 
can set a = 77 = 5 = 0, then it can be shown that i? = 
and that the resulting spacetime is fiat. Also, if one triad 
leg Ca is chosen according to Eq. (|5.ip . it is not possible 
to apply a boost to render the null vector geodesic, or 
equivalently to set the coefficient /3 to zero. An exam- 
ple which illustrates this fact is in the asymptotic region 
of a radiating spacetime, where our choice of c° would 
point along a cylindrical radius; meanwhile, the outgo- 
ing null geodesies define a radial direction, and it is clear 
that these two directions are not orthogonal. Rather, we 
would need to locally choose some other null direction to 
define Z". 

We now ask whether it is possible to find a coordinate t 
whose gradient is timelike and orthogonal to Cq, i.e. that 
t^ac" = 0. The first step is to define a timelikc unit vector 
Ta as 



1 
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{la + ria)- 



(5.10) 



From the normalization conditions (j3.3p it is straight- 
forward to verify that TaT" = —1. We would like to 
determine if Ta is hypersurface orthogonal, so that it can 
be written as Ta = —Qt_a- This is possible if and only if 
Ta is twist-free, T\^aDbTc\ = 0. In 3D, this is equivalent 
to the vanishing of the scalar 

W = e-'Wclb = ^i-p + j-C + 2rj + e). (5.11) 



For a general l" and n° this will not be true, but we can 
choose a boost A that will transform 77 such that W = 0. 
By Eq. (|5.9|) . we see we must choose 



1 



?7 



(/3-7 + C-e)- 



(5.12) 



We have so far fixed our triad, and selected the harmonic 
coordinate "0 and the coordinate t whose gradient lies 
parallel to Ta- Let us call the third coordinate s. On the 
coordinate basis (t, s, ip) the assumptions thus far imply 
that in all generality the we can express the covariant 
components of the triad as 



h = {—hi hs, ~h^) , 
c, = (0, 0, v^). 



{~nt, ~hs, h^) 

(5.13) 



The factor g in the definition of Ta \s g = {It + nt)/V2. 
The metric on the coordinate basis is constructed using 
Eq. p.5|) . To see if any further metric functions can be 
set to zero, consider a coordinate transformation that 
leaves the coordinates t and tp unchanged but chooses a 
new coordinate s', such that s = f{t, s' , ip). We find that 
the metric can be expressed in the same form except with 
the functions h, rit, hs, transformed as 



hsf,t , 



h' 



n 



hsf.s' , 



(5.14) 



It is thus always possible to choose a gauge in which 
/t,^, — 0. Dropping the primes, the resulting metric on 
the coordinate basis is 



hij = I -hs {k - nt) 



-2ltnt -hs {k - nt) 
2hl 

I 



(5.15) 



For the rest of this section we make this coordinate 
choice. The covariant components of the triad vectors 
are 



r 



1 



V2gh, 



■{hs, It, 0), n' 



1 



V2gh, 



■{hs, -nt, 0) . 



= (0, 0, /r72). 



(5.16) 



The choice oi %p as a, coordinate is an unfamiliar 
one, and to help build some intuition we present the 
Minkowski metric, triad, and rotation coefficients in this 
coordinate system in Appendix O The rotation coef- 
ficients in general axisymmetric spacetimes can be ex- 
pressed in terms of the functions entering Eqs. (|5.13p 
and (|5.16p . and are listed in Appendix [P] 

The expression in (|Dip for the coefficient e. 



(5.17) 



2^/2 

can be integrated using the Bianchi identity (|5.2p . 

^ {\nR)^^ 
^ 4V2 



(5.18) 
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Combining these equations shows that the metric func- 



tions obey [\n R{hgg)^ 
provides 



0, which after integration 



ih,g)^R^ g{t,s). 



(5.19) 



There is stih some residual coordinate freedom in 
Eq. ()5.15p in that we can apply a coordinate transfor- 
mation to the s and t coordinates without changing the 
form of the metric. In particular by using the coordi- 
nate transformation s — /2(^, and using a restricted 
version of Eq. (|5.14p it is possible to choose a gauge in 
which g{t,s) = 1 so that we have [hgofR ^ 1. We will 
not necessarily make this specialization in the rest of the 
text. 



1. Field equations adapted to the ip coordinate choice 

In this subsection we specialize the field equa- 
tions p.l5p ~ p.23p to the case where we use as a coor- 
dinate and where Ca = \J2/R'ij,\a- Recall that this choice 
implies that 7 = — e, V,i = V',2 = and V,3 = vW^- 
With this specialization, we re-order the general field 
equations p.l5|) - p.23p augmented by the commutation 
relations (|5.3p . One of the field equations is redundant 
with one of the commutation relations, while the remain- 
ing 11 equations can be split into a subset of four equa- 
tions that contain directional derivatives in the l"^ and 
directions only, 

e.i = +Pa + 2/3(5, 
Ci --e,2 + 2aC, 
L_i = —0_2 + at — 59, 

5 1 = a,2 - i?/2 + 2a(5 + /3C + (5.20) 

and a group of seven equations that fix the directional 
derivatives of certain rotation coefficients in the c° direc- 
tion, 

e^3 = —(5,1 — t.i + a_2 + 2a(5 + ab + 9b + 2e^, 
P,3^0.i+0{a-9) + 2P{e-7-i), 

= — oiri + ae + /3((5 — b) — rj9 — 6e, 
(5,3 ^ ri^2 - oe( + S{r] + e) + C9 + V- - '-e, 
C,3 = ^2 -Si + 20] + b^ + 2Ce, 
9^3 = 2e,i -I3b- rj9 + e{-e), 

^3 = -2e,2 + C^ + t(??-e). (5.21) 

We showed that in the coordinate basis (i, s, i}}) the met- 
ric can be written in the form (|5.15p . With the choice of 
1°- and n° in (|5.16p . all the equations (|5.20p contain only 
derivatives with respect t and s, and effectively constitute 
a set of constraint equations that have to be satisfied for 
every constant ^ surface. 

As can be seen from the above set of equations, choos- 
ing -0 as a coordinate does not greatly simplify the field 



equations. For this coordinate and triad choice the ma- 
jor simplifications occur in the Bianchi identities (|5.2p . 
the form of the metric ()5.15p , and the simple form of the 
corresponding Weyl scalars. An additional advantage of 
this coordinate and triad choice that will be discussed in 
the next section is the easy identification of the axis. 



2. Axis conditions as A — >■ 

On the axis, which for the three metric is denoted by 
the boundary conditions ij: — > —00 or A — > 0, we now 
explore the conditions on the triad quantities required 
for the elementary fiatness condition to hold. 

The first step is to observe that working in a coordinate 
system where "0 is a coordinate makes it easy to prove the 
equivalence of the two forms of the axis conditions, Kjj = 
1 in Eq. (|2.28p and the coordinate invariant expression 
in Eq. (j2.29p . Assuming the metric hij can be written 
in the form ()5.15p , the metric on the space orthogonal to 
the axis 1^2 is merely ds^ = g^^ dijP' + Xd(fP, where the 
4D metric component g.^^ = 2/(Ai?). The elementary 
fiatness condition ()2.28p now reads 



lim — r- 

J —o 



1, 



(5.22) 



where use has been made of the fact that A is not a func- 
tion of 0. Applying rHopital's rule and differentiating 
above and below the line with respect to ip, the elemen- 
tary fiatness condition becomes 



lim ■ 

A->0 



-—^ = lim 

/g4>ip 



(5.23) 



or equivalently R —J- 2e~^''f'. By definition, R = 2'0,a'0'°j 
showing that the covariant expression (|2.30p and thus 
Eq. (|2.29p are equivalent to the elementary flatness condi- 
tion. Note that the elementary fiatness condition, in con- 
junction with the condition found when examining the 
rotation coefficient e, Eq. (|5.19p . implies that the deter- 
minant of metric on the subspace normal to the axis also 
remains finite as we approach the axis. To see this explic- 
itly, observe that det[ft,^j] — — 2/i^p^, where i,j G {s, t}. 
By the condition found in Eq. (|5.19p in the gauge where 
g{s,t) = 1 we have det[ft,^j] = —2/R. The determinant 
associated with the corresponding part of the four metric 
becomes det[gj-] = —2/{Re^'^), which by the elementary 
flatness condition approaches the value —1 on the axis as 
expected. 

Symmetry dictates that a null vector on the axis re- 
mains on the axis when it is sent out to infinity or 
toward the origin. Thus on the axis and are 
geodesic, provided they are chosen to lie along the in- 
going and out-going directions. In terms of the NP 
scalars (|4.12p . this translates into k — ^er^ l\/2 — > 0, 
and V = Ce3'>/^/2 0. 
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In Appendix [Dl explicit formulas for the expansions 
of the metric quantities about the axis are given and 
discussed. The special case of the static Schwarzschild 
black hole is examined in Sec. I VII CI where the scaling 
of the solution, the 3-curvature R and all the rotation 
coefficients are explicitly computed. 



B. Geodesic null coordinates 

We now examine the equations in a coordinate system 
adapted to asymptotic null infinity, where the concept of 
emitted radiation is well defined. Akin to the standard 
methods used in the NP formulation (see e.g. [2^, HI, 
losj). this coordinate system and triad choice is tied to 
the tangent vectors of null geodesies. We begin with a 
family of null hypersurfaces in S, and we label these by a 
coordinate u, so that h^^u^iUj = 0. We then choose the 
covariant representation of one null triad vector to be the 
gradient of the coordinate u, setting li = —u^i. Since U is 
the gradient of a coordinate, it has vanishing curl. The 
intrinsic derivative of li on the triad basis is 



'a|6 





a 





■n 



(5.24) 



The fact that la,\b is symmetric immediately sets /3 = a = 
and r] = —e. Note that /? = a = implies that la is 
geodesic and affinely parametrized on S, by Eqs. (|3.7|) . 
Also, recall that null geodesies are conformally invariant, 
and we can verify that here Eqs. (|4.12p imply that if /3 = 
then K = in A^. Thus if 1°" is the generator of a geodesic 
null congruence on S, the corresponding null congruence 
in the physical manifold is also geodesic. The above con- 
ditions on the rotation coefficients further imply that the 
field equation p.22p is trivially satisfied. If we choose 
as another coordinate the affine parameter p along the 
geodesic that l'^ is tangent to, we have l^f-^i = /_i = dpf. 
Lastly, we label our third coordinate Expressing the 
null vectors on the (w, p, x) coordinate system, we have 



= (-1,0,0), r = (0,1,0). 



(5.25) 



The normalization conditions (j3.3p allow us to restrict 
some of the components of the remaining triad vectors, 
giving n" = 1 and = 0. Using the expression for the 
metric in terms of the triad vectors p.5|) . we can see that 
f^uu ^ f^ux and = -1 follows. Three more met- 
ric functions fully determine hij. We parametrize these 
remai ning metric components following the convention 
of [9^ Il00l | so that the contravariant form of the metric 
becomes 







-1 



-1 2vi4 
V2e' 





V2e^ 



'2V3 



(5.26) 



The covariant form of the metric on S is then given by 

/ ~2vi -1 V2e^^ \ 
h,j =(-10 I . (5.27) 



V2e" 







This metric holds for any null foliation of the manifold S, 
where constant u surfaces denote the null hypersurfaces, 
the affine parameter p serves as a coordinate along a par- 
ticular geodesic and the coordinate X; usually associated 
with an angular coordinate, labels the geodesies within 
the hypersurface. 

We further need to fix the triad legs n* and c*. 
One such choice that satisfies the normalization condi- 
tion p.3p and gives the correct form of the metric (|5.26p 
is 

= (l,-t-i,0), c' = iO,V2,e-^'). (5.28) 

The corresponding covariant vectors are 

n, = {-vi,-l,V2e^'), c, = (0,0,6"^). (5.29) 

On this triad, the directional derivatives applied to a 
function / are 



= f,p , J. 2 = J.u - Vlf^p , /,3 = t'2/,p 



(5.30) 



If the chosen coordinates (u, p, x) are to be valid, 
they must satisfy the commutation relations given in 
Eq. p.29|) . Applying the commutation relations to each 
successive coordinate provides a simple way of relating 
the rotation coefficients to derivatives of the metric func- 
tions of Eq. (j5.27p . The commutators acting on x yield 
the coefficients 



7 = -e, 



= -""3,1, 



V3.2- 



(5.31) 



Applying the commutation relations to u reiterates that 
a = /3 = and rj = — e. Finally, applying the commuta- 
tion relations to p fix 



e = ^ {V20- V2,l) , 

C = V2V3,2 + V2,2 + Wl,3- 



(5.32) 
(5.33) 
(5.34) 



1. Field equations adapted to the geodesic null coordinate 
choice 

When working with geodesic null coordinates, where 
J = 1] = and a = /? = 0, the field equations p.lSp - 
p.23|) can be expressed in simplified form in terms of the 
five remaining rotation coefficients as 



(5.35) 
(5.36) 
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'5,1 


= e - Y"</'4, 


(5.37) 


f ,2 


^ f ,3 


— OU — c — UL , 

2 


yo.oo j 




+ 0.2 


= -OS, 


(5.39) 


Ci 


+ e,2 


= 01 , 


(5.40) 




- t,2 


= /,('- - S) +(f>2, 


(5.41) 


■3 4 


-2e,2 


= Ce-2ie. 


(5.42) 



Ci + 5. 



One of the equations is trivially solved and has been 
omitted. The remaining equations have been reordered, 
and some are linear combinations of the original set. 
These combinations are (|07)) = ([3^/2 - ((XT5]) /2 
— p.l9p describing the derivative (5_i; the combina- 
tion ([Og|) = ((X^ /2 - ((XT5|) /2. yielding an expression 
for the combination 9^2 — e,3; and finally the combina- 
tion (|5.40p = p.23p — (|3.16p to obtain an expression for 
Ci + e.2- The remaining equations are simplified ana- 
logues of their counterparts in Eqs. p.l5p - p.23|) . 

The reordering makes apparent the fact that a hier- 
archy exists in the reduced system of equations, which 
in turn makes it possible to formally integrate the field 
equations in a systematic way. Suppose we begin on 
a null hypersurface of constant u on which the direc- 
tional derivatives of the function -0 are given, so that (pi , 
i ~ ... 5 arc known. Equations (j5.35p and (|5.31|) can 
be integrated with respect to p to obtain the rotation co- 
efficient 9, and subsequently the metric function v^. In a 
similar fashion Eqs. (|5.36p and (|5.34p yield e and V2, and 
subsequently (|5.37p and (|5.32p give 5 and vi . The metric 
functions W3, V2 and vi are thus determined within the 
null hypersurface up to boundary terms. The require- 
ment that D'^ifj = determines t = ^3^2 using Eq. (|3.30p . 
Thus the manner in which the metric function W3 changes 
away from the initial null hypersurface is known. Equa- 
tion (|5.38p then serves as a consistency condition which 
restricts some of the six integration constants that arise 
while integrating Eqs. (|5.35p - (|5.37p . The other integra- 
tion constants are determined by boundary conditions 
that will be discussed more fully in Sec. IV CI Equa- 
tion (j5.40p implicitly determines W2,2- Equation (|5.4ip . 
in conjunction with the condition D'^ip = provides 
an evolution equation of ip. The remaining two equa- 
tions, (|5.39p and (|5.42p . are eliminant relations that are 
trivially solved when the metric functions are substituted 
into the field equations. 

The hierarchy of field equations that arise when they 
are expressed on a coordinate system adapted to a null 
hypersurface has been extensively studied in the four di- 
mensional context. It is known, for instance, that the 
equivalent equations on M. are for mally integrable on a 
constant u surface [s^ [9^ 

[io1[1Q1- The asymptotic be- 
havior of the metric in geodesic null coordinates, and the 
associated boundary conditions are further discussed in 
Sec. IV CI where the relationship to the Bondi formalism 
is explored. In Appendix |F] we give an explicit example 
of how the field equations are systematically integrated 
in an asymptotic region far from a gravitating system, 



although there I"' is affinely parametrized with respect to 
M. 

The gauge and triad choice discussed in this section 
has the advantage of eliminating four of the nine rotation 
coefficients. This reduction in complexity makes appar- 
ent a hierarchy in the field equations that hints at the 
possibility of finding an analytic solution to the prob- 
lem. In Sec. IVI Al we carry out an example calculation 
in which the field equations (|5.35p - (|5.42p arc systemati- 
cally solved in a special case. It should be noted however 
that this simplicity comes at a cost. Unlike the case 
where the triad was adapted to the coordinate tp, and 
the Ricci tensor only had one non-zero component, the 
Ricci tensor on this triad is constructed from the three 
independent quantities ip^a- It has two degenerate eigen- 
vectors with zero eigenvalues, and a single normalized 
eigenvector with nonzero eigenvalue 2ip^a'4''°'- 

The analysis performed in this section assumes that 
li is affinely parametrized in S. If we adjust the pa- 
rameter along each geodesic, p — > p'{u,p,x)i this re- 
sults in hup' 7^—1, but otherwise preserves the form of 
the metric. A physically motivated alternative to affine 
parametrization in S is to boost /' so that it is affinely 
parametrized in the physical 4D spacetime, and then to 
use the affine parameter r as the second coordinate in- 
stead of the parameter p. The field equations that result 
form this choice of parametrization are detailed in Ap- 
pendix |e1 

Some gauge freedom remains when P is affine in S, 
and p is used as a coordinate. Shifting the origin of 
the affine parameter along each geodesic separately, p' = 
p -j- f{u,x)^ transforms the metric function of (|5.27l) ac- 
cording to 



V2 



V3 ■ 



(5.43) 



Relabeling the individual geodesies within a spatial slice, 
x' = fl(u,x), transforms the metric functions of (|5.27p 
to [lOl 



V2 



e ^V2 -e ■ 

9,x ' 



9,n 

9,x 



9,x ' 
(5.44) 



Finally, it is also possible to relabel the null hypersur- 
faces, setting u' = h[u), p' ~ p/h^^- The metric compo- 
nents transform as 



^, v'^ = V3. (5.45) 



C. Asymptotic flatness and the peeling property 

We will complete our discussion of useful coordinate 
systems on S by discussing the asymptotic limit of the 
metric far from an isolated, gravitating system. We will 
consider only spacctimes that are asymptotically fiat and 
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therefore admit the peehng property |82H8a ]. According 
to the peehng property, the Weyl scalars expressed on an 
afhnely parametrized out-going nuU geodesic tetrad ad- 
mit a power series expansion at future nuU infinity (de- 
noted 1+ ) of the form 



(") 



(5.46) 



n=0 



where are constant along an out-going geodesic, 

i.e. ^'i"''(M,x)- The work of Bondi, van der Burg and 
Metzne r 19911 . as weU as Tamburino and Winicour's ap- 
proach [98[, indicate that the the metric functions also 
admit a power series expansion if expressed in terms of 
geodesic null coordinates. 

Appendix |F] details a triad-based derivation of the 
asymptotic series expansions of the metric functions, in 
the restricted context of axisymmetric spacctimcs. In 
this derivation, the "Bondi news function" is identified 
with the derivative of the dominant coefficient in the ex- 
pansion of the shear of the out-going null tetrad leg. The 
calculation is performed assuming that the out-going null 
geodesic is affinely parametrized in M, and the corre- 
sponding field equations given in Appendix |E] are used. 

The results obtained in Appendix iFl for the asymptotic 
expansion of the metric can be summarized as follows. In 
terms of affinely parametrized null coordinates, the 4D 
line element can be expressed as 



2du dr 



5^'^W2 du dx 

(5.47) 

where, according to Eqs. ((FTT|) . ((FT5|) . ((FTsI) . 

and (|F2ip . the metric functions admit the following 
asymptotic expansion as the affine parameter t — > oo: 



r(0)2 



(0)~ 



3t 



+ 0(r-2), 



Wl 



r(0)2 



(i-x^) 



W2 



[(1 



1 
2^2 

(0)1 



6r2 



(0) 



(0) 



(1 - x2)t2 



^2 ^3 



(5.48) 

Note that the coordinate x is chosen here to be x = 
cos 9, and 9 is the usual polar angle. The axis occurs as 
X^±l- 

The free functions that enter into the metric arc 
a^^^ {u, x) and the dominant terms associated with the 
Weyl scalars ^f\u, x): * G {0,1,2}. The dominant 



terms of ^"3°^ and '^^l'' arc fixed by these free functions 
through Eqs. (|F25|) and ((F26)) of equivalently 



^4 " ,uu ■ 



(5.49) 



The field equations determine the evolution of '^f^ (ti, x) , 
i € {0,1,2} from one null hypersurface to another via 
Eqs. (|F22|) . (|F27| . and (|F28)) . As can be observed from 
(|5.49p . the free function a'^u' carries the gravitational 
wave content of the spacetime and is often referred to 
as the "Bondi news function." 

A solution that settles down to a Schwarzschild black 
hole in its final state requires that in the limit u — >■ 00 
the scalars behave as 

{^-0, *2, ^-3, *4, <T<"' , a|°^} ^ {0, 0, -A/, 0, 0, 0, 0} , 

(5.50) 

where the constant M is the mass of the final black hole. 
For u < 00, 4'-°^ , i G {0,1,2} are then determined by 
these final conditions, provided that cr*-"-' (u, x) is given, 
using the evolution equations (|F22|) . (|F27|) . and (|F28p . 
For easy reference these equations are repeated here : 



= 3a(°)M") + (l-X^) 



(0) 



= 2a(")*f + 
(0) [(1-X^)-<°\J 



(0) 



2 a 



2^1 -X^ 



(5.51) 
(5.52) 



2 



CT^"V("\„„, (5.53) 



We now examine how the metric and Weyl scalars be- 
have on the axis in the limit of large distance from the 
isolated source. As noted in Sec. IIIFl and explored fur- 
ther in Appendix IdJ the metric functions have a power 
series expansion in A = e'^'^ near the axis of symmetry. In 
addition, these expansions are such that the metric func- 
tions vanish sufficiently quickly in the approach to the 
axis so that there are no "kinks" at the axis [9^. Note 
that from (|5.48p , the coefficient of the du dx term in the 
metric, namely e'^^W2, is only regular on the axis if both 
a^^^ and '^^^ vanish on the axis, and respectively scale 
like 



rW = (l-x')a(°) 



(0) 



(5.54) 



where ct*^°^ and ^J"'' need not vanish on the axis. Sub- 
stituting these scalings into the evolution equations for 
the dominant expansion terms for the Weyl scalars, 
Eqs. (|5.51l) - (j5.53p and (|5.49p . shows that on the axis 



vl^o = *! = *3 = *4 = , 



(5.55) 



indicating that the spacetime is Type D on the axis, and 
there is no radiation to infinity along the axis. This is 
to be expected, since spin-2 transverse radiation cannot 
propagate along the axis and still obey axisymmetry. The 
only nonzero Weyl scalar is 4*2, and the dominant coef- 
ficient can depend only on u. 



^-2 



-M{u)t~ 



(5.56) 
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The metric functions in the near-axis, large r limit are 
2M{uy 



9uu 



1 



(5.57) 



with the g^y. term becoming singular at the poles simply 
due to our coordinate choice. Changing from x to the 
coordinate 9 gives gee = + 0[t) while fixing gue = 
on the axis. Asymptotically, the only dynamics present 
are the variation of the multipole moments with changing 
M, where M(u) clearly gives a monopole mass moment. 

The results given thus far are for a metric whose r 
coordinate coincides with the afhne parameter of the 
geodesic null vector l^^ on the physical manifold M.. In 
order to convert to afhne geodesic null coordinates on 
the manifold S, and so read off the asymptotic behav- 
ior of the metric (|5.27p . we need to consider the effect 
of the transformation between (u,p, x) and (w, r, x) co- 
ordinates, where p ~ p{u,T,x)- Expanding dp in ()5.27|) 
in terms of du, dr and dx and equating the result with 
(|E3[) yields the following relationship between the metric 
functions and the derivatives of the afhne parameter p, 



V2 = e^'^'W2 +P,x^~ 



V3 = W3, 

vi = e '^wi ~p.u- 



(5.58) 



Integrating the first equation of (|5.58p with respect to r 
yields 



2^|^ -(0)^2 ^ ^(0)2^ ^^jj^. 



P =(1 - X ) - cr''"r" +cr^' '-T , g 

+ 0(r-i). 



(5.59) 



Inverting the series (|5.59p to obtain an explicit expression 
for T in terms of p is complicated by the logarithmic term. 
The leading order expression can however easily be found 
and is 



VI. SOLUTIONS TO THE TWIST-FREE 
AXISYMMETRIC VACUUM FIELD EQUATIONS 

In this section we will characterize the properties of 
known twist-free solutions to the axisymmetric vacuum 
field equations within the framework that was established 
in the previous sections. The aim is to identify existing 
solutions and to catalog the assumptions made in finding 
them. With the exception of the Schwarzschild solution, 
none of the existing asymptotically fiat solutions have 
physical significance. The hope is that this characteri- 
zation will help establish the necessary properties a new 
dynamical solution, such as the head-on collision, must 
posses. 

A number of insights that can be gleaned by relating 
the four dimensional physical quantities to the three di- 
mensional rotation coefficients are discussed in Sec. IIIII 
This section thus relies heavily on Sec. lIVl and in particu- 
lar Eqs. (|4.12p and Eqs. (|4.15p . which give the 4D NP spin 
coefficients and associated Weyl tensor in terms of the 3D 
rotation coefficients discussed in Sec. IIIII Wherever pos- 
sible we will also express the properties of the known so- 
lutions in terms of the two geometrically motivated triad 
and coordinate choices of Sec. IV Al and Sec. IV Bl 

We begin the discussion of analytic solutions with an 
example of the systematic solution of the field equa- 
tions mentioned in Sec. IV Bl We consider the special 
case where the spacetime admits a coordinate choice in 
which we can simultaneously choose as a coordinate 
with Ca = y^2/Rip,a and find a null coordinate u such 
that the geodesic null vector la = —u a is orthogonal to 
Ca- This example has the benefit that it draws on our 
general results for both coordinate choices discussed in 
SeclVZland Sec.lVB] 

Having found one solution, we then place it in context 
with known solutions using a classification scheme based 
on the optical properties of the geodesic null congruence 
that 1°' is tangent to. It should be noted that the scope 
of many of the known solutions discussed in this section 
often extends beyond the restricted arena of twist-free 
axisymmetry, but we will restrict our discussion to this 
realm. 



By working out the series expansions of Eq. (|5.58p in 
terms of t and then substituting Eq. (|5.60p into the re- 
sult, the asymptotic behavior of the metric (|5.27p can 
be found to be 

-1=^(1 -xY/^(l + 2a(°)) 

+ {ipY/%l-X^f'^^f'^ +0{\np) 



2/3 



2a(") 



^\'\0{p-'''\np) 



e"^ =(3p)2/3(l - x')i/3 + o(p-i/3 inp) . (5.61) 



A. Special case: spacetimes that admit the 
coordinate choice (it, p, ijj) 

Any three metric can be expressed on an affinely 
parametrized geodesic null coordinate basis (u, p, x) as 
in Eq. (|5.27p . In this section we will consider the special 
case where the third coordinate x = V'- Making the triad 
choice defined in Eqs. (|5.25p and (|5.29p we thus require 
that ijj^i be spacelike and orthogonal to li. From the re- 
sults in Sec. lVBI on the geodesic nuU coordinate choice we 
have that a = /3 = 0,e = — 77 = — 7 and that the simpli- 
fied field equations presented in Eqs. (|5.35p - (|5.42p hold. 
As discussed in Sec. IV Al the choice of as a coordinate 
naturally sets ip^i = -0^2 = and the metric function 
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g-2i;3 — Rj2. Furthermore the only non-zero curvature 
scalar is <j)Q ^ R which also simplifies Eqs. (|5.35p - (|5.42p . 
We now proceed to solve this set of field equations. 

First, we note that Eq. (|5.35p . 6,p = 9^ can be solved 
by setting = — [p + f{u, ip)]^-^. We use the coordinate 
freedom discussed in Sec. lVBI to relabel the origin of the 
affine parameter p by an arbitrary function of /(it, -0), to 
give 



-P 



.1) 



Performing one more integration using the commutation 
relation (|5.3ip . v^^p = p~^, allows us to obtain the metric 
function v^; equivalently, the scalar curvature R = 2e^^"^ 
is 

R = cxM)p-''. (6.2) 



The next field equation (|5.36p . e^p = indicates that 
e = e(u, V') only. Once more a commutation relation can 
be integrated to obtain the metric function vi. In this 
case Eq. (|5.33p . {v2p),p = — 2ep implies that 



V2 = C2{u,1p)p 



ep. 



(6.3) 



Before proceeding, let us use the fact that ip has 
been chosen as a coordinate and examine the simplified 
Bianchi identities (j5.2p . The first equation is trivially 
satisfied, while the third equation (lni?)^3 — 4e places re- 
strictions on the integration constants already obtained. 
Writing out the directional derivative in terms of coordi- 
nate derivatives and substituting in the solutions for i?, 
V3, e and V2 we have 

2e--^2^==p-i+2c2(u,V')p-'=0. (6.4) 
2^ci(u, V) 

This expression must vanish for all powers of p, which 
implies that e = C2 = ci_^ = 0. A consequence of this 
result is that V2 = 0, and in addition 7y = 7 = e = 0, 
and ci ~ ci{u) is a function of u only. The final Bianchi 
identity gives the coefficient t. 



Cl. 



2ci(u) 



vip 



(6.5) 



Substituting the results obtained thus far into the 
third field equation, ()5.37p we obtain d^p = — i?/2, and 
thus the rotation coefficient S = ci(u)/(2p) -I- cs{u,tp). 
The integration constant cs{u,'>p) is fixed to the value 
cs = ci.„/(2ci) by evaluating the field equation (|5.38p . 
The commutation equation (j5.32p . vi p = —5, yields an 
expression for the final metric function vi, 



Vl 



ci{u) 



Inp ^ 



Cl.u 

2ci 



■p -f- €3(7/, V') • (6.6) 



Since V2 ~ 0, the final commutation equation (|5.34p sets 
C = C3^,0-\/ci/2p^-^ but the field equation (|5.40p implies 
that = 0. Thus we have that C3 = 03(7/) is a function 



of u only and C = 0. All metric functions now depend on 
the variables p and u only. 

It is useful to observe that we still have the freedom 
to relabel the null hypersurfaces of constant u as dis- 
cussed in Eq. (|5.44p . If we transform to a new set of 
coordinates {u', p',ip) such that u' ~ J du^J ci(u)/2 and 
p' = p-\/2/ci(w), the metric function v'^ expressed on 
the new coordinate basis can be written in the form 
v'l = — Inp' + c'^{u'). It is always possible to choose a 
coordinate u that labels the null hypersurfaces in a man- 
ner such that ci(u) = 2. For the rest of the section we 
make this choice (omitting the primes). 

The final field equation (|5.4ip reduces to = — i(t— 5). 
Upon substituting in t = ~vip~^, 6 ^ p^^ and vi — 
— Inp + C3 we find that C3.„ ~ 0, so that if we define 
hiA = C3,A is merely a constant. 

In summary, when i/; j is orthogonal to a geodesic null 
vector li = u^i the metric functions are 

vi^ln(^^^, V2=0, V3^\np, R ^ 2p-\ (6.7) 

and the rotation coefficients take on the values 



7/ = 7 = C = 0, 



= -- i = -^ 6^- 



P 



P 



P 



(6.8) 



Having successfully solved the 3D field equations in 
this special case, let us examine some of the implications 
the solution has for the 4D spacetime associated with 
the original axisymmetric problem. The 4D NP spin co- 
efficients for the solution found in this section are easily 
obtain from Eqs. ()4.12p 

k — i) — I3^i~0, f — — 7T = — 



V2p- 

These expressions for the spin coefficients show the corre- 
sponding 4D null congruence is also geodesic and affinely 
parametrized. By writing down the Weyl scalars using 
Eq. dSJl) 



"^2 



p^ 



p3^ 



V2p' 



P 



(6.10) 



we observe that /'^ is a principle null direction. A gen- 
eral classification scheme using the spin coefficients will 
be discussed more fully in Sec. IVICI For now it is use- 
ful to observe that the fact that p = a and that l'^ is a 
geodesic principle null vector indicates that this space- 
time is a cylindrical-type Newman- Tamburino solution 
10J|. These spacetimes do not depend on any free func- 
tions of w. In fact, the metric of Eq. (j6.7p corresponds 
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to the particular case of a cylindrical-type Newman- 
Tamburino spacetime with one of the two arbitrary con- 
stants that parametrize these solution set to zero [2g . 

The axis conditions offer no additional constraints to 
this solution. As we approach the axis, we have e'^'^ — >■ 
0. In order for the axis to be free of singularities, the 
elementary flatness condition, Eq. p.29p must hold. In 
this particular case, we would have h^^'4',i'^,j = P"^) ^'Hd 
so elementary flatness would require 



lim e'^'f' =p, 



as we approach the axis. However, we then have that 
p — >■ as we approach the axis, and we can see from 
the Weyl scalars ()6.10p that the spacetime is singular as 
p — > 0. We can conclude that this solution possesses a 
curvature singularity along the axis. 



B. Spacetimes with special optical properties 

In the next subsection we review known, special solu- 
tions to the axisymmetric, vacuum field equations, clas- 
sifying them according to their optical properties. The 
classification will be made according to the properties of 
the null congruence that the tetrad vector is tangent 
to in , and by extension the congruence that the triad 
vector P is tangent to in S. One of the benefits of the NP 
formalism is that the spin coefficients are directly related 
to the optical properties of a given spacetime. By seeking 
solutions with specified optical properties, many simpli- 
fications become possible, and the assumptions made are 
physically transparent. As detailed in Sec IIVI by choos- 
ing to work in twist-free axisymmetry we have at least 
halved the complexity of the problem of solving the 4D 
field equations. The 4D spin coefficients computed from 
the 3D rotation coefficients are all real, which has imme- 
diate implications for the null congruence they describe 
in 4D, and we will discuss these implications here. 

Consider a general axisymmetric spacetime whose ax- 
ial KV is twist-free and explore the behavior of the con- 
gruence of null curves that is tangent to in A^. The 
expansion and twist of this congruence are described by 
the real and imaginary parts of the spin coefficient p, 
respectively, and constitute the first two optical scalars. 
From Eqs. (|4J2)) we know that 



3.12) 



and is manifestly real, and this shows that a spacetime 
with a twist-free, axial KV admits a twist-free null con- 
gruence (c.f. [261 llOSl p. Furthermore, a twist-free con- 
gruence is hypersurface orthogonal, and so the fact that 
p is real implies that lu is proportional to the gradient of 
some potential function u. As in Sec. IV Bl in this case 
is geodesic, since I'^lfj.-i, = w^u-^u — (w'''u,^)^^,/2 ~ 0, 
and we have 



Twist-free Axial KV 

k = ^\p\ = 



PND geodesic , 



4-0 = 



, PND not geodesic 



*o 7^ 



a = 



3 j^-) Type II, RT AxisymNT 



p = 



Bilge 



Physical 
spacetimes 



Kundt 



FIG. 2: Classification of spacetimes possessing a twist-free, 
axial Killing vector. The abbreviations used in this figure can 
be interpreted as follows; Killing vector (KV), Principle null 
direction (PND), Robinson- Trautman (RT), and Newman- 
Tamburino (NT). 



The final optical scalar that characterizes the geometri- 
cal properties of the null congruence is the shear, which 
measures the distortion of the congruence and is given 



by 



3.14) 



K = e-'^/3/y2 = 0. 



(6.13) 



An interesting property that arises from restricting the 
discussion to twist- free axisymmetric spacetimes is that, 
if the spacetime further has a null direction along which 
the derivative of ip vanishes, the associated congruence 
has p ^ a. The vanishing of any of the other directional 
derivatives of "0 gives analogous reductions to the 3D ro- 
tation coefficients, as can be seen by studying Eqs. (|4.12p . 

A number of solutions to the field equations have been 
found which admit a geodesic, hypersurface orthogonal 
null congruence (where = 0). We will discuss these 
spacetimes and their relation to the form of the field 
equations developed in this paper in the next subsection. 
Our focus will be on asymptotically fiat solutions, which 
can represent isolated systems, and we will reserve our 
discussion of stationary, axisymmetric spacetimes until 
Section IVIII Figure [2] gives a summary of the solutions 
we will be considering, along with the reductions and 
assumptions employed to yield the known results. 



C. Principle null geodesic congruences 

We have showed that in any twist-free, axisymmet- 
ric spacetime there exists a geodesic, hypersurface- 
orthogonal null congruence. If in addition, the tangent 
to this congruence is also assumed to be a principle null 
direction, so that \l/o = 0, the field equations can be 
solved and the exact metric expressions arc known. To 
see this, it is easiest to work in a triad where is affinely 
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parametrized with respect to the physical manifold A4. 
In this case, we have a = —2ip^i on 5, e = on A^, and 
Eqs. (|4TT5l) and (|3?20| become 

«'o = = V^2^+Vaii, (6-15) 
e,i = {9 + ijjsf (6.16) 

Note that the directional derivative in these equations 
can be interpreted derivative with respect to the 

affine parameter r in and expressed as /.i = /.r- For 
Eq. (|6.15p the two solutions are 

V-,! - or V.i = 7 ■ (6.17) 

T + Cl[U, X) 

Using Eq. (|6.15p . Eq. (|6.16p can be rewritten as 

(6* + 1 = (^ + V^.i) • The two solutions to this equa- 
tion are 

6' + Vi=0 or 6l + Vi = r. (6.18) 

T + C2{U,X) 

By substituting the solutions (|6.18p and (|6.17l) into the 
4D spin coefficients, Eqs. (|4.12p . several distinct condi- 
tions on the optical scalars p and a can be identified. 
The expansion free case, where 9 — ijj^i ~ 0, implies 
that p = a = Q and is known as the Kundt solution. 
On the other hand if tp^i = 0, or if ^? = we have 

that p = ±(T ^ which characterizes a cylindrical-type 
Newman- Tamburino spacetime. In the case that V',! 7^ 
and t/' i + 7^ the transformation r — > t' + f(u, x) can 
always be used to set ci = —02- The optical scalars thus 
become p = —t/{t'^ — c\) and a = ci/(t^ — cf). This 
case can be split into two distinct scenarios: If ci = 0, 
(7 = and the spacetime can be classified as a Robinson- 
Trautman spacetime. If on the other hand ci 7^ 0, then 
a is nonzero and the result would again be a Newman- 
Tamburino spacetime, but of spherical type. However, in 
the case of axisymmetry, we can solve the field equations 
explicitly for a nonzero ci , by integrating the hierarchy of 
field equations and matching powers (and transcenden- 
tal functions) of t at each step; the resulting spacetime 
has vanishing curvature and so is actually flat. This con- 
forms to the known fact that the (non-trivial) spherical 
Newman- Tamburino solution can have at most only a 
single ignorable coordinate, namely the parameter label- 
ing the null hypersurfaces u [26l H 06] : in other words, the 
spherical-type solutions are incompatible with axisym- 
metry. 

In the subsequent subsections we examine the prop- 
erties of the each of the spacetimes mentioned here in 
greater depth. 

1. Newman-Tamburino spacetimes 

Newman- Tamburino spacetimes are characterized by 
the properties 

K = *o = , =0 and a^O (6.19) 



The met ric for these solutions can be found explicitly 
jl03lll04{ and except for special cases, the spacetimes are 
of the generic Petrov Type I. The Newman-Tamburino 
solutions are divided into two classes, "spherical type" 
and "cylindrical type" solutions. The spherical type is 
the more general, requiring 7^ &a* . The cylindrical 
type requires p^ = aa*. O nly t he cylindrical type solu- 
tions admit a spatial KV [l06l |. and so are the case of 
interest for our study. Since all of the 4D spin coeffi- 
cients are real in twist-free axisymmetry with our tetrad 
choice, these solutions require that the more restrictive 
condition a = ±p holds, or equivalently in terms of the 
3D quantities, 

either ip,i = or 6* = -ips- (6.20) 

In the case with ip^i = 0, the proper circumference 
of the orbits of the axial KV is unchanging along the 
geodesic null congruence. Therefore, these solutions rep- 
resent a spacetime that expands in the direction of the 
congruence. Note that the congruence is not simply 
frozen at a constant parameter r, since the expansion 
p is nonzero for these solutions. 

The Newman-Tamburino solutions do not correspond 
to spacetimes of physical interest. In addition, the metric 
functions have a simple polynomial dependence on the 
coordinate u, which shows that the dynamics of these 
spacetimes are very simple. Si nce the pr operties of these 
solutions are well understood |l03l Il04l | and the general 
derivation of the metric functions is lengthy, we do dis- 
cuss these spacetimes further. Instead, recall that the 
solution found Sec. IVI Al ls a special case of the cylindri- 
cal type Newman-Tamburino solutions, where in addition 
to V',! = we assumed that 'ijj^2 = 0. The solution found 
in Sec. IVI Al is parametrized by one constant A, while 
the genera l cylindr i cal-t ype metric contains two arbitrary 
constants ^ [Tol, [To!] . 

2. Robins on- Trautman spacetimes 

The second class of solutions where the congruence is 
geodesic, principle null, shear-free, (k = it = ^0 = 0) and 
expanding ( 8 ^ 0) is known as the Robinson- Trautman 
spacetimes jlOTl llOSf . The solutions to the field equa- 
tions in such a case can be reduced to a single nonlinear 
partial di ffere ntial equation and have been well studied, 
see e.g. (109{ and the references therein. While these 
equations have been used to study radiating sources in 
an exact, strong field setting, they do not represent phys- 
ical systems, such as a stage of head-on collision of black 
holes. 

In terms of the 3D rotation coefficients, the conditions 
for the Robinson Trautman solutions are 

/3 = 0, 6* = -2^/;,i 7^ 0. (6.21) 

Note that since these vacuum spacetimes admit a shear- 
free geodesic null congruence, the Goldberg-Sachs theo- 
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rem (25l. [26l . [TToj states that they are algebraically special, 

so that *o = *i = 0. 

We can verify this directly from cm the 3D equations. 
To do so, we use an affine parametrization with respect 
to the physical manifold A^, as discussed in Appendix lEl 
This sets a = —2ip^i. Substituting 6 = — 2?/'.i into the 
field Eq. (|6.16p immediately gives — 0. Showing that 
VPi = can also be made to vanish requires more finesse. 
For this, we use Eq. (|4.15p with f3 — to obtain an 
expression for 

*i = ^[^,i(V',3 + ?7)+V',i3]- (6.22) 

In Eq. (|6.22[) the commutation relation p.29|) has been 
used to interchange to order of differentiation on ip. The 
field equation (|3.18p . specialized to the case where 6 = 
a = — 2'0_i and /3 = 0, can be written as 

r/4 = -2[V'4(Vj,3 + 2r/) + V.is]- (6.23) 

Observe that with the simplifications so far (k ^ a = 
e = *o = 0), the two 4D Bianchi identities [11, Ull that 
govern only the directional derivatives of '^q and ^'i can 
be expressed in our notation as 

J = 4/5^1 3 = (2/3 + 47^)^1 (6.24) 

The corresponding expressions in terms of the 3D 
quantities associated with the triad discussed in Ap- 
pendix [EJ where e = —7 = — r/ — 2?/; 3, is 

ln(e'"^^'i),i = 0, ln(e4'^l'i),3 = -5tj. (6.25) 

If wc apply the directional derivative l^di to the sec- 
ond of Eqs. (|6.25p . use the commutation relations p.29p 
to switch the directional derivatives, and repeatedly use 
Eqs. (|6.25p . we obtain the equation 77^1 = —2r]ip^i. Substi- 
tution of this into Eq. (j6.23p implies that V'.i ("0,3 + v) = 
—^.13, and thus ^1 = 0. 

To complete the discussion of the Goldberg-Sachs the- 
orem for twist-free axisymmetric spacetimes, note that 
substituting the condition ^0 = = into the 4D 
Bianchi identities immediately implies that k = a ~ 0, 
and thus the existence of a geodesic shear-fee null con- 
gruence. 

3. Expansion-free spacetimes 

Wc now consider the last case in our class of spacetimes 
that admit a geodesic principle null congruence, namely 
spacetimes that are both shear-free and expansion-free 
(k = a = p ~ ^0 ^ 0). These metrics are Kundt s olu- 
tions and have been extensively studied [2^ lllll . Ill2| . In 
terms of the 3D quantities, the Kundt metrics have the 
properties j3 = 9 = ip^i = 0. Since ip^i = 0, setting a = 
implies that the geodesies can be affinely parametrized 



in both M and S simultaneously. The fact that = 
also greatly simplifies the 3D curvature scalars, setting 
05 = 04 = 03 = 0. Observe that choosing li = —u.i and 
setting 9 = in addition to a = (3 = and jy = — e im- 
plies that D^u = 0, so that u is a harmonic coordinate. 
This can be seen by taking the trace of Eq. (|5.24p . 

In the Kundt metrics, many of the rotation coefficients 
and metric functions arc independent of the affinc param- 
eter p, which simplifies the calculations. The solution of 
the field equations in this case provides another simple 
and illustrative example of the integration of the 3D field 
equations. We now proceed to solve the hierarchy of field 
equations (|5.35[) - (|5.42p in Sec. IV Bl in conjunction with 
the commutation relations applied to coordinates. 

In addition to V'.i = 0, which implies that ijj is indepen- 
dent of the affine parameter, Eqs. (|5.36p and (|5.39p give 
e^i = i.i = 0. Further, the commutation relation (|5.3ip 
yields 1)3 _i =0. Thus ^3, ip, e and l are functions of only 
two coordinates u and x- Note that it is always possible 
to use the coordinate freedom to define a new coordinate 
x' = g{iJ-,x) such that = in the new coordinate sys- 
tem, by choosing — e^-'. Since l — V3^2 by Eq. (|5.3ip . 
it is possible to set 

V3 = L = 0. (6.26) 

The next metric function, V2, can be found using the 
commutation relation (|5.32p . V2, p = — 2e. Since e = 
e(it, x) only we can integrate the equation to yield 
1)2 ~ —2ep + ci(u,x)- The ability to shift the origin of 
the affine parameter by a function of u and x by defining 
a new parameter p' ^ p + g{u, x) allows us to set ci = 0, 
and thus 

V2 = ~2e{u,x)p- (6.27) 

In obtaining the expressions for W3 and V2 we have used 
up most of the coordinate freedom with respect to the x 
and p coordinates, except for transformations of the type 
x' = X + /i (fi) and p' ^ p + f2 (u) which do not spoil any 
of the simplifications so far. 

In order to learn more about the function e consider the 
harmonic condition on tjj, Eq. (|3.30p . With the reductions 
employed thus far, Eq. p.30p reduces to ip^^^ = 2eip,3- 
Furthermore the field equation (|5.38p simplifies to the 
expression 6^3 = + %p\. Taking the smrr and the differ- 
ence of these two equations we obtain 

(e + V',3),3 = (e + V',3)^ (e - V',3),3 = (e - V',3)^. 

(6.28) 

For functions of the form / = .f{u,x) the fact that 773 = 0, 
implies that the directional derivative /_3 becomes the 
coordinate derivative /.;^, for these functions it is further 

true that /,2 = 

Before solving Eq. (p?^ . note that Eq. ((07)) can be 
reduced to S,i = — indicating that i5.i is indepen- 
dent of p and allowing us to integrate the equation to 
find an explicit expression for S, 

S^ie^-^%)p~wsiu,x), (6.29) 
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where wg is an integration constant. In addition, S has to 
satisfy the difference of Eqs. (|5.42|) and (|5.40p . 6 3 + e.2 = 
— 2?/'^2V',3- Substituting in Eq. (|6.29|) and evaluating, we 
obtain the following constraint on ws'. 



(6.30) 



Integrating the commutation relation (|5.32l) . 5 = — 
the metric function vi can found to have the form 



Vl 



ip%)p'^ +ws{u,x)p + wi{u,x) 



3.31) 



where Wi is a new integration constant. To see what 
additional constraints are to be imposed on the integra- 
tion constants by the field equations, we express the only 
remaining coefficient C in terms of the metric functions 
using Eq. ([53i|) with U3 = 0, 
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Substituting this expression for C into Eq. (|5.4ip and ex- 
pressing the result using coordinate derivatives yields the 
constraint 



(6.33) 



This is the final condition that has to be satisfied. 

All that remains now is to explicitly integrate 
Eq. (|6.28p . There are three possible cases: 



e = 



1 



and 
and 

and 



^,3 



0, 



(6.34) 



1 



C2{u) ' 

(6.35) 

1 



(6.36) 



where the functions Ci(u) are arbitrary functions of u 
only. For the remainder of the discussion we shall con- 
centrate on the most generic case, Eqs. (|6.36p . Taking 
the difference of the Eqs. (|6.36p and integrating one more 
time gives an expression for ip, 



The Ricci scalar R 



In 



X - C3(lt) 



X-C2(lt) 



C4{u). 



3.37) 



2ijj'^ and the coefficient e are 



R = 



C3 - C2 



ix - C2){X - C3) 

C3 + C2- 2x 
ix- C2){x- C3) 



(6.38) 
(6.39) 



It is straightforward to verify that the solution given 
here matches that presente d in the original derivation 
of Kramer and Neugebauer [l05l | , who also use the con- 
formal 2-t-l decomposition. The solution provided by 



Hoenselaers using the triad method [2J| fails f or the twist- 
free case, but Mcintosh and Arianrhod jll2l | show (after 
making some corrections) that it matches the above form. 

We now investigate the condition of elementary flat- 
ness to see if it provides any additional constraints on 
the free functions. The axis is located at those points 
where 



„2V 



,2C4 



X-C3 
X- C2 



0, 



(6.40) 



which occurs when x — ^ C3. The condition for elementary 
flatness, Eq IJ^TI^ or equivalently e'^'^i? 2, can be 
expressed as 



o4c4 



1™ \' 

A^O 4(C3 - C2)' 



= 1 , 



(6.41) 



The requirement that the axis be free of conical singu- 
larities does thus provide a further constraint on the free 
functions, relating C4 to C3 — C2. 

Let us further explore the properties of the spacetime 
by computing the Weyl scalars. Since the Kundt space- 
time is geodesic and shear free, it is algebraically special, 
with ^0 = ^1 = 0. The other Weyl scalars can be ob- 
tained in terms of the triad variables by making use of 
Eqs. Km . 



(6.42) 



^ — (?/.^x„ + eV;,„ + 3^,„V,x) ' (6-43) 

*4 = e^'^ (i^,„„ + 3(^,„)2 - - CV'.x) , (6-44) 

This shows that the metric is in general of Petrov Type 
II. In the case where 02 = 03, we can sec from our solution 
in Eq. (|6.37p that = C4(u), and so using Eq. (|6.42p we 
have that ^2 = 0. In this case the metric is of Type III. 
In case 2 mentioned in Eq. (j6.35p . e = —ip^ ^^^'^ implies 
^2 = and that the spacetime is also of Type III. Fi- 
nally, we consider the special case 1 of Eq. (|6.34p where 
e = -0^^ ~ 0, for which the vanishing of the Weyl scalars 
^2 = *I'3 = implies that the metric is Type N. This 
completes the full classiflcation of axisymmetric space- 
times that admit a geodesic principle null congruence. 



D. Other axisymmetric, twist-free spacetimes 

Given the complete classification of spacetimes with 
special optical properties, it is clear that dynamical 
spacetimes of physical interest are not represented in this 
class of solutions. We must necessarily consider space- 
times whose geodesic, hypersurface-orthogonal congru- 
ence is not a principle null congruence, and 4*0 7^ 0. 

Known exact solutions that do not admit a geodesic 
principle null congruence include the twist-free solutions 
with the restriction (p/ir);^/^ = which have been found 
by Bilge |ll3| . However, Bilge and Giirses also showed 
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that this class of spacGtime s, th ough generally of Type I, 
is not asymptotically flat jll4l |. The class of twist-free 
spacetimes that obey {p/a)-fj.l^ = includes the vacuum 
Generalized Kerr-Schild (GKS) metrics, when they have 
twist-free congruences. GKS metrics are of the form 

g^.L' = .9a"' + HiJ„ . (6.45) 

where g^^, and g^^, are both solutions to the vacuum field 
equations; is a geodesic null vector with respect to 
both metrics (it forms the twis t-free congruence); and H 
is some function on spacetime jll3l |. Gergely and Perjes 
showed that the GKS spacetimes which admit twist-free 
congrue nces are the homogeneous, anisotropic Kasner so- 
lutions |115| (although two of the three constants that 
classify the Kasner spacetime must be set equal in order 
for there to be a single rotational symmetry and so an 
axial KV). 

Another method to find solutions with one spatial KV, 
which we mention for completeness, is to use a different 
triad choice than those discussed here. The idea, as de- 
veloped by Perjes in the study of stationary spacetimes 
[4l|, is to orient one of the triad legs along an eigenray. 
The eigenray is defined such that if its spatial tangent 
vector is geodesic in the spatial 3D space, it is the pro- 
jection of a null ray in the full 4D space. The triad for- 
malism of Perjes adapts naturally to the case of a spatial 
KV rather than a timelike KV, in which case the eigenray 
is timelike and the triad is chosen such that 

^,i=V',2, ^,3 = 0. (6.46) 

If the eigenray vector is, in fact, assumed to be geodesic, 
then solutions to the field equations can be found . In 
stationary spacetimes these were found in [4l|, Ill6| , and 
in the case of a s pacet ime with a spatial KV they were 
found by Lukacs jll7l |. When the geodesic eigenray is 
shearing, the solutions are Kasner (as in the case of the 
twist-free GKS spacetimes), and when it is not shearing 
they are Type D, and so very restricted. 

Thus, future studies which aim to extract physical in- 
formation about isolated dynamical, axisymmetric space- 
times will have to focus on general spacetimes, where 
none of the principle null directions are geodesic, and 
which do not fall within Bilge's class of metrics. 

VII. STATIONARY AXISYMMETRIC VACUUM 
SPACETIMES 

As discussed in Sec. HI the SAV field equations have 
been completely solved, in the sense that techniques ex- 
ist that can generate any SAV solution. In this section, 
we briefly discuss the case of SAV spacetimes in the con- 
text of the conformal 3D metric and Ernst equation of 
Sec. Ill D[ We then specifically focus on the twist-free case, 
and catalog the simplifications to the triad formalism 
discussed in Sec. IIIII when static spacetimes are consid- 
ered. We recast the static metric in the form discussed 



in Sec. IV Al with tp chosen as a coordinate. Finally, as 
an illustrative example, we explore the properties of the 
Schwarzschild metric re-expressed using ip as a coordi- 
nate, and we derive the scaling of the rotation coefficients 
as the axis is approached. We compare the results to the 
general expansions derived in the time-dependent case in 
Appendix |DJ 

A. SAV spacetimes with twist 

A stationary, axisymmetric spacetime possesses two 
KVs. One is spacelike, with closed orbits, which we de- 
note ^'^ as in previous sections. The other is timelike, 
and we call it rj^. The ignorable coordinate associated 
with 77^ we denote t, and we will later work in a gauge 
were rf = 5]. The KV 77^ places additional restrictions 
on the spacetimes considered up until now. 

Carter showed that for axisymmetric spacetimes which 
arc asymptotically flat, the two KVs commute [o^l- In 
addition, in vacuum, the pair of KVs are surface forming, 
which allows us to use coordinates t and cj) in order to 
separate the metric into a pair of two dimensional blocks. 
Since any two metric is conformally flat, we can further 
choose coordinates p and z such that these coordinates 
are isotropic on their block [89| . and express the metric 
in the Weyl canonical form, 

ds^ =e-2^ [e2'= {dp" + dz^) + p^dcP"] - e^^idt - Adc^f , 

(7.1) 

where the functions U, A, k arc functions of p and z only. 

This form of the metric is extensively used for explor- 
ing the SAV field equations. Associated with the metric 
functions is an Ernst potential, which is often used in 
constructing solutions to the field equations. It should 
be noted that the Ernst potential usually employed in 
the discussion of SAV spacetimes in the literature is 
associated with the timelike KV 77'^ , and is not the Ernst 
potential £ associated with the axial KV introduced 
in Eq. ((2?23)) . In the SAV context, = e^^ + iip where 
e^^ = —77^77^ and (^^^ = efj^apV'^V^''^ ■ In order to make 
a direct connection with the notation used in this paper, 
we redefine the metric functions and cast the metric of 
Eq. (fTTjl in the form 

ds'- =e^^ {dcp - Bdtf + [e^^ {dp^ + dz^) - pHt^] . 

(7.2) 

with e^"^ = A = £,^£,fj, as before. The metric functions of 
the two forms are related by 

Using the metric (|7.2p . the conformal 3-metric is given 
by 

/i,, =diag[V,e2^e^^]. (7.4) 
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It turns out that for the hue element (j7.2p . the Ernst 
equation for f = X + iuj and field equations for "0,7, 
and B can be solved in an identical manner to the more 
usual SAV case, where U, k, and A are sought and when 
f(^) = / + 'i'P- As derived previously, the field equations 
reduce to a single nonlinear equation for £, namely the 
first equation in Eqs. (|2.23l) . Written out in the coor- 
dinates associated with metric (|7.2p . the equation for £ 
becomes 



10 - 



£ + £* 



(7.5) 



Here we have defined as the usual fiat space Laplace 
operator, in cylindrical coordinates. When working with 
SAV spacetimes it is often useful to introduce the com- 
plex coordinate C = {p -'r iz) / to express the equations 
more compactly on the complex plane. In these coordi- 
nates = {dp — idz)/V^- Once £ is known the metric 
functions for 7 and B can be found by making use of the 
line integrals 



7.C 



V2p{£)^c{n. 



C 



{£ + £*f 
2p{£^£*),C 
{£ + £*Y 



(7.6) 
(7.7) 



B. Twist-free SAV spacetimes 

We now specialize to the case of static axisymmetric 
spacetimes, which arc twist-free. The KV is hyper- 
surface orthogonal, and thus lo ^ A = B = {). The Ernst 
equation (j7.5|) reduces to the 3D cylindrical Laplacian 
applied to ijj, 



VV = /o"^(/oV',p),p + V,.. =0. 



The two metric functions 7 and '0 are related to the met- 
ric functions in Wcyl canonical metric (|7.1|) by 



Tp — \n p — U , 



7 = lnp + fc - 2f/. 



(7.9) 



Since Inp is a homogeneous solution to the cylindrical 
Laplacian, we see that if U obeys the cylindrical Laplace 
equation, then does also, and vice versa. All homo- 
geneous solutions to the cylindrical Laplace equation are 
known (for example in terms of Lcgendrc polynomials). 
The only difficulty is in specifying boundary conditions 
whose corresponding solution gives a spacetime with the 
desired physical interpretation. We will examine this is- 
sue using specific examples. 

In the context of the line element (j7.ip . a single 
Schwarzschild black hole is generated by a line charge of 
length 2M placed on the axis. With this as the boundary 
condition, the Laplace equation for U can be solved; in 
Weyl coordinates, the solution is 



^2U ^ r-+ + r_ - 2M 
r+ + r_ + 2M ' 



(7.10) 
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FIG. 3: Orthogonal coordinates -0 (solid lines) and s (dashed 
lines) plotted for a single black hole of mass M = 1 in canon- 
ical Weyl coordinates (p, 2). 



2'Y 2 

e ' = p 



(r+ + r_ + 2Mf 
4(r+ + r_ - 2M)r+r. 



(7.11) 



where r'i^^ p^ + {z± M)^ and e^"^ = p'^e"'^^ . The hmc- 
tion ip for the Schwarzschild metric in (p, z) coordinates 
is plotted in Fig. [3l The black hole lies on the axis be- 
tween z /M = ± 1 , where the equipotential lines of "ip meet 
the axis at almost right angles. Further from the axis the 
surfaces of constant ■0 rapidly approach surfaces of con- 
stant cylindrical radius. For the single black hole, as we 
will discuss in the next section, the appropriate axis con- 
ditions guaranteeing elementary flatness arc satisfied at 
p = outside the line charge. 

Multiple static black holes solutions can be found by 
placing multiple line charges along the axis |ll8l | . Mathe- 
matically this corresponds to the superposition of multi- 
ple Ui potentials given in Eq. (|7.10p . centered at positions 
Zi and with masses Mi. The corresponding potential 
is then constructed using ip = Inp — X^I^i^*- Given 
the Tp potential of a superimposed set of black holes, we 
can always construct a solution for the metric function 7 
satisfying Eq. (j7.6|) . which for static spacetimes becomes 



(7.12) 



Despite the fact that a solution can be found, it is not 
possible to find a solution for which elementary flat- 
ness holds along every connected component of the axis. 
From Eq. (|7.ip we see that if fc 7^ along a compo- 
nent of the axis, then elementary fiatness does not hold 
there. Similarly, comparison of the line element (|7.2p 
with Minkowski space written with ip as a coordinate [see 
Eq. (jCip in Appendix [C] . shows that elementary fiatness 
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requires j = 2'ip — \np along the axis. The conical singu- 
larities that result when the elementary flatness condition 
is not met are interpreted as massless strings or struts 
which hold the black holes apart, keeping them station- 
ary. The fact that these singularities always appear in 
static, axisymmetric black hole solutions is in line with 
our intuition that black holes should attract each other, 
and can never remain stationary at a fixed separation. 
Approximate models based on using the tension on the 
strut to evolve bi naries in a head-on collision scenario can 
be found e.g. in jll9l - ll2l| . In the time-dependent case 
the harmonic equation governing ?/; once again suggests 
that a generalized superposition principle could hold, at 
least on the initial time slice. 

In the time-dependent case, however the counterpart 
of equations (|7.12p defining the gradient of the potential 
7 do not exist. Instead the second order elliptic equations 
for 7 associated with the initial value problem must be 
solved. This allows the freedom to impose the boundary 
conditions guaranteeing elementary flatness. We thus ex- 
pect that in dynamical spacetimes the conical singular- 
ities can be removed. We have yet to fully explore the 
time-dependent equations in the framework provided by 
this paper. 



Static Axisymmetric spacetimes with ip chosen 
as a coordinate 



In order to explore SAV spacetimes within the frame- 
work of more general axisymmetric spacetimes, we trans- 
form to a coordinate and triad system adapted to the 
scalar field ip. To do this, consider the coordinate sys- 
tem (t, s, if)), where t remains the ignorable coordinate 
associated with the timelike KV, -tj} \s & potential that 
obeys Eq. (|7.8|) and the coordinate s is orthogonal to il^ 
and defined by 



dpS 



dzS = -pdpip. 



(7.13) 



The integrability condition for s is guaranteed by the 
vanishing of the Laplace equation (|7.8p for tp. Changing 
coordinates from (t, p, z) to (t, s, V') yields a metric in 
the form 



h,j = diag[-p2, 2/5, 2/R] 



(7.14) 



where the functions S and R are the normalization fac- 
tors defined by S" = 2s;iS'' and R — 2ip.iip''^. Note that 
from the definitions in Eq. (|7.13|) and the metric (|7.4p . 
it follows that S = p^R. This is consistent with the 
result (|5.19p obtained by integrating the equations for 
the coefficient e, in the special case where g{t,s) = 1. 
This means that once R and p{s, ip) are known the en- 
tire metric in ("0, s) coordinates is known. In terms 
of the Weyl canonical coordinates (p, z), we have that 
R = 2(?/'^p + ?/»^2)e~^'*', and in addition from the static 
field equations for 7, Eq. (|7.12[) . we also have the iden- 
tity R' = 4e-^^^l + jl)/p\ 



Let us now take a closer look at the rotation coeffi- 
cients associated with the metric (|7.14[) . When we con- 
sider the more general metric in Eq. (|5.15p . make the 
triad choice (|5.16p . and substitute both that It = rit and 
that the metric functions are independent of time into 
Eqs. (|DT|) . we find 



(7.15) 



a = (5, 7 = -e, ^ = -C, ?7 = 0, 



Furthermore note that for the choice of triad vec- 
tors (|5.16p . the directional derivatives obey the relation 
/2 = — /,! for any time-independent function /. 

Noting that = l/(p^i?) and that g = p in the static 
case, the remaining independent rotation coefficients can 
be expressed in terms of only the functions _R, p and their 
derivatives as follows. 



R 



4:V2R 

pR.s 

4Vi?' 



/3 = - 



Rp,. 



a 



V2p 

Rp.s 



(7.16) 



To further illustrate the implications of choosing ip as 
a coordinate, we now turn to the concrete case of the 
Schwarzschild metric. Computationally it is useful to 
introduce prolate spheroidal coordinates (x, y) related to 
the Weyl coordinates by the transformation 

p^ ^ M^{x^ -l){l^y^), z^Mxy. (7.17) 

In spheroidal coordinates, D'^ip = is equivalent to re- 
quiring that 

[{x^ - 1)9,^] + dy [(1 - y^)dytP] = 0. (7.18) 

In terms of the {x, y) coordinates, the norms of the 
Schwarzschild azimuthal and timelike KVs are 



{x + \f{l-y\ 



x + l 



Note that the upper and lower segments of the axis are 
identified by the coordinate values y ~ \ and y — 
respectively. The event horizon of the black hole is indi- 
cated by X = 1. By direct substitution it is easy to verify 
that both V and U satisfy Eq. ([71^ . 

It is further possible to verify that the potential s = 
(x — \)y has a gradient orthogonal to the gradient of ip 
and obeys Eq. (|7.13p . It can thus be used as the second 
spatial coordinate. The lines of constant ip and s coordi- 
nates are plotted in the (p, z) plane in Fig. |31 Since in 
Weyl coordinates the [p, z) plane is conformally flat, the 
fact that the curves intersect orthogonally in Fig. [3] indi- 
cates that their gradients are orthogonal to each other. 
The strong warping influence of the black hole at p = 0, 
— 1 < z < 1 on constant ip surfaces in these isotropic 
coordinates is clearly visible. This behavior can be as- 
cribed in large part to the coordinates, which compress 
the black hole horizon onto the axis. As e''^ — )■ away 
from the hole, contours of constant -0 approach the axis. 
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but on the black hole tp acts as an angular coordinate, 
changing as the surface is traversed. 

The 3D Ricci curvature scalar R for the Schwarzschild 
metric is 



R 



2{x + 2j/2 - 1) 



2e 



(a- + l)= 



(7.19) 



and obeys the axis condition (|2.30p . Also note that 
Eq. (j7.19p is written is the same form as the more general 
series expansion of R about the axis given in Eq. (|D8|) . In 
Schwarzschild case, the series truncates after the first or- 
der. To facilitate compact notation later on, let us define 
a function Rq that is finite on the axis by 



Rq 



R 2y2 



1 



(7.20) 



The nontrivial rotation coefficients for the 
Schwarzschild metric are 



Se^^jx - 1) 
2Rl^\x + l)e 



3e'^V^^v/(a: + 1)2 ~ 
V2i?;^/'(.T + 1)11/2 



a2lp 



e''*(3a;-7)- 1 



^3/2 V 2(2; + 1)6 2(a; + l)V ' 



^/2i^o/'V^^(x + 1)5/2 



(7.21) 



From the above expressions it is clear that on the axis, 
e2'> = {y = ±1), we have that (3 = 6 = while e 
and a diverge. The fact that /? = is an indication 
that the null vector 1°^ is geodesic on the axis, as is ex- 
pected from symmetry. The divergence of a is an indi- 
cation that this geodesic has been poorly parametrized. 
A better choice would be to boost the null vector so that 
it is affincly parametrized. Note that in the boosted 
frame the vector (l" + n'')/V2 = {Al" + A~^n'')/V2 
is no longer hypersurface orthogonal. Choosing A so 
that the component of the vector along the t coor- 
dinate direction corresponds to that in the Kinnersley 
frame |l22| results in an affine parametrization on the 
axis. This boost transformation is achieved by setting 
A = p{x + l)/[\/2{x — 1)]. The transformation of the ro- 
tation coefficients into the boosted, affinely parametrized 
frame are given in Eq. (|B2[) . and is straightforward to 
compute. 

In this frame, the Weyl scalars computed from 
Eq. (1121) are 



*2 



3e2'A 



32^(33; + 1) 

' Ro{x + l)^ ^ 2i?o(a; + l)6 



y 



4Ro{x + 1)5 
1 



2V2i?o(a 



3e'''(a;- l)y 
V2Roix + 1)5 



*4 



3e2V'(:r-l)2 
" Roix + iy 



(7.22) 



On the axis, the only nonzero Weyl scalar is 
^2 = ~(a; -t- 1)"'^, where x + 1 can be associated with 
the standard Schwarzschild radius along the axis. As 
we move off the axis the other Weyl scalars take on 
non-zero values. This is expected because our choice of 
triad is adapted to the gradient of tp rather than being 
a geodesic null triad. In fact, anywhere off the axis in 
the Schwarzschild spacetime, a triad adapted to the tp 
coordinate can never have its null vector be geodesic, 
as can be verified using Eq. (j7.2ip . This provides an 
explicit example of the general arguments regarding the 
boost transforms in (|5.9p . 

When the Schwarzschild metric is recast into a form 
with tI) as a coordinate, many of the equations appear 
to be unwieldy and offer little additional insight, but 
they do allow for an explicit verification that a frame 
exists in which the near-axis scaling of the rotation co- 
efficients and metric functions computed in Appendix iDl 
hold. The main motivation for choosing a triad adapted 
to the ip coordinate is that the results obtained in this 
analysis of the SAV case generalize to dynamic space- 
times. The Weyl scalars and rotation coefficients should 
give some indication of the expected behavior of their 
counterparts in dynamical spacetimes, near the axis and 
near black holes. It should be noted that many of the 
features that make the standard SAV analysis elegant 
are due to the availability of isotropic (p, z) coordinates, 
and the ability to linearly superimpose multiple solutions 
in the static case. These properties do not generalize to 
time-dependent spacetimes. One feature that is common 
to both the SAV and the time-dependent analysis is that 
■0 is a harmonic function that plays a crucial role in de- 
termining the configuration of the spacetime. It is hoped 
that a generalized superposition principle by which a new 
solutions for tp can be formed by the "sum" of two or 
more existing solutions can be found. Such a superposi- 
tion of solutions, although straightforward to achieve on 
an initial value slice, will subsequently be complicated by 
the fact that the potential tp influences the evolution of 
the metric functions that determine its own evolution in a 
possibly nonlinear way, making a "sum" of two solutions 
nontrivial. 



VIII. CONCLUSIONS 

In this paper wc have reviewed the reduction of the 
Einstein field equations in the case where the spacetime 
admits an axial Killing vector. The problem of find- 
ing solutions to the field equations in 4D then reduces 
to a problem of finding solutions to the field equations 
in 3D with an additional scalar field source term. We 
specialized to the case of vacuum spacetimes, and then 
to twist-frec spacetimes, where the field equations be- 
come especially simple, but are still capable of describing 
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spacetimes of physical relevance. Of particular interest 
is the case of a head-on collision of non-spinning black 
holes. In order to recast the equations into a form that 
seems especially amenable to investigation and intuition, 
we have presented a triad-based formulation of the equa- 
tions due to Hoenselaers. We have expanded upon the 
original work of Hoenselaers, linking this formalism to the 
Newman-Pcnrosc formalism and discussing two triad and 
coordinate choices that help to simplify the equations. 
We have also reviewed the known twist-free axisymmet- 
ric solutions which are not necessarily captured by the 
SAV equations, classifying them according to their opti- 
cal properties, which correspond to certain assumptions 
simplifications in the field equations. 

We have also introduced and explored the use of a har- 
monic coordinate ^ on the 3 manifold S. Recall that S is 
conformally related to the manifold S of orbits of the KV. 
The function ip corresponds to a scalar field source for the 
Einstein field equations on 5, and its use as a coordinate 
simplifies the curvature on S considerably. There is a 
natural expansion about the axis of symmetry in terms 
of A = e^"^' , which provides inner boundary conditions 
for the field equations on S, and we have provided these 
series expansions and their connection to the elementary 
flatness condition. We have also revisited the case of 
static, axisymmetric spacetimes using -0 as a coordinate, 
in order to concretely illustrate the near-axis behavior of 
the triad and curvature quantities. Meanwhile, the as- 
sumption of asymptotic flatness provides the usual outer 
boundary behavior for the metric functions and curvature 
quantities, in terms of a Bondi expansion in geodesic null 
coordinates far from the sources. In order to develop a 
uniflcd notation, we have given explicit expressions for 
the Bondi expansion in the triad formalism, eventually 
arriving at series expansions for the metric hij on S. 

We intend this work to serve as a comprehensive and 
usable reference for the challenging goal of arriving at 
new solutions to the fleld equations. In future work, 
we hope to investigate initial data schemes, both null 
and spatial, which may help to explore the problem of 
the head-on collision using analytic techniques. Another 
spacetime of interest which is dcscribablc by the equa- 
tions reviewed here is the case of critical collapse of grav- 
itational waves in axisymmetry; this case likely has ad- 
ditional symmetries, as is seen in the numeric studies of 
critical collapse [H, [t^], which can help guide further 
investigation of this solution. 

It is also possible that some extension of the solution 
generating techniques which allowed for the solution of 
the SAV equations can be modified and carried over to 
the case of the axisymmetric equations. The exploration 
of how the techniques developed to solve the Ernst equa- 
tion might be extended to dynamic spacetimes will be 
the subject of future study. 
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Appendix A: Dimensional reduction of the 4D field 
equations 

In this appendix we review some of the results pertain- 
ing to the curvature of the three dimensional manifold S 
whose induced metric /i^^ is related to the metric on the 
four dimensional manifold A4 by Eq. (j2.3|) . In this sub- 
section ^'^ is not necessarily a KV, but merely assumed 
to be timelike, i.e. ^^^^ = A > 0. Just as in Sec. [Hi 
Eq. (j2.6p , the covariant derivative operator Da is defined 
by the full contraction of the 4D derivative operator with 
the projector /i" = (5" — A"-'^^^^". For convenience, the 
definition of Da , Eq. (|2.6|) is repeated below, 

Dafp^ = hi^Xp^iV ^%p) . (Al) 

1. Generalized Gauss-Codazzi equations for a 
timelike projected manifold 

As in the case of the 3-1-1 split in numerical relativ- 
ity [131, the Gauss-Codazzi equations describe the rela- 
tionship between the 3D and 4D curvature tensors asso- 
ciated with the metrics /i^^ and respectively. The 
Gauss equation can be derived by considering the 3D 
Ricci identity, which defines the contraction of the 3D 
Riemann tensor with an arbitrary covector Va on 5, 

RaMsV" ^'^D^'yD&]yc.- (A2) 

The derivation proceeds by writing out the derivative op- 
erators on the right hand side in terms of the 4D quan- 
tities g^u^^i^L and using Eqs. (|2.3p and (|Aip . As an 
intermediate step we define the quantity Kpa to be 

Kpa = -A-i/2/,M/j^v^e„ (A3) 

and expand the double covariant derivative operator ap- 
plied to Vy as 

DaDpVy =K^hPV^, [KhlV^Vr) 

+Kpah:,\-^'H''^.Vr. (A4) 
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To derive the second line of Eq. (|A4[) wc repeatedly use 
the definition of K^p given in Eq. (|A3p . the fact that 
hi^h'^ = and the identity Vr^ ^ . 

Substituting Eq. (|A4[) into Eq. (jA2p results in an ex- 
pression relating the Riemann tensor on S to the Rie- 
mann tensor on A4, 

+ 2K[s^^hlX-'^^^^V,Vr. (A5) 

This result holds for any vector with norm A. The 
top line of Eq. (jA5P resembles the usual Gauss equation 
often encountered in a 3+1 split of spacetime if the ten- 
sor Kai3 is identified with the extrinsic curvature of the 
embedded hypersurface (there is a relative sign change 
in front of the terms containing quadratic products in 
the tensor Kpa that results from the fact that we are 
considering a timelike rather than spacelike 3 manifold) . 
The tensor Kap defined in Eq. (|A3p can be identified 
with the extrinsic curvature of a hypersurface embedded 
in A4 only if is hypersurface orthogonal. The second 
line of Eq. (|A5[) contains a term with the prefactor K[s^] . 
In general if ^'^ has twist, this term is non- vanishing and 
must be retained. 

For all vectors V" whose Lie derivative with respect to 
vanishes, C^Va — 0, wc can simplify the second term 
in Eq. ^ using hlX~^/^^''VaVr = RfSaV^ to yield 

+ 2K[s^]Kp^V^. (A6) 

Since the vector is arbitrary provided C^Va =0 and 
V^^'^ = 0, it can be dropped from Eq. (|A6|) to give an 
expression for the curvature on the 3 manifold in terms of 
projected quantities. An important consequence of the 
3D Riemann tensor so obtained is that in order for it 
to have the correct symmetries the tensor KjSa must be 
antisymmetric, Kpa = ^[/3q] • The condition K(^aip) = 
is the same as requiring that the projection of the 4D 
Killing equation hold. 

A concise way of expressing the generalized Gauss 
equation Eq. (|A6p in terms of the vector ^'^ in the case 
where -fS^(a^) = is 

Ro^p^s =h'^h''phPh''sR 

+ ^hl^^h^h^,^ (V^e(.) (Vp)^.) . (A7) 

Note that the since the curvature tensor Rap-yS is defined 
on a 3D manifold which has a vanishing Weyl tensor, 
Rap-yS can be constructed solely from the Ricci tensor 
Rap- As a result only the contracted Gauss equation 

RpS =h'phs {Rvcr — ^^^S,'^^'^ R^vpa) 

+ \^'''^J^''p]^,^ (Vm^i.) (Vp)^.) , (as) 

needs to be considered. 



The components of the 4D curvature tensor where one 
index has been projected onto n^, the unit normal in the 
direction (and not to be confused with the null triad 
or tetrad vectors used elsewhere in this text), are related 
to quantities defined on the 3 manifold S via the Codazzi 
equations. These equations can be derived by applying 
the 4D Ricci identity to the unit vector — A"^/^^'^ 
and projecting the result onto the 3D manifold S 

KKh1n''R,,,p„ = 2KKh1VypV,}n^ . (A9) 

When expanding the right hand side of Eq. (jA9p in terms 
of 3D quantities, it is useful to observe that Kpa can be 
expressed as the gradient of the unit vector n'' which has 
been twice contracted with the the projection operator. 
Expanding this relation and using the fact that n^n^ — 1 
yields the identity 

V^?i^ = -Kyf, + Uf^a^, (AlO) 

where ap = ^np is a measure of how the unit vec- 
tor n'^ is changing when parallel propagated. Note that 
haO-tJ. = since a^n'* = 0. Substituting Eq. (jA10|) into 
Eq. (jA9p yields a generalized Codazzi equation 

h'^hPhgn''R^^p^ = DgKa^ - D^K^s + 2aaK[^iY 

(All) 

The last term once again vanishes in the case where <^'' 
is hypersurface orthogonal, but has to be retained if we 
consider vectors with twist. Contracting Eq. (|A11[) on 
the indices a, 7 with the metric h"^ yields the contracted 
Codazzi relation 

h^sn^R,, - DsK - + 2r%5], (A12) 

where K = K^. 

There is one more nonzero contraction of the 4D Rie- 
mann tensor with the projector and the unit vector 
71^, which is computed by contracting the second and 
third indices of the Riemann tensor with and the re- 
maining indices with the projection operators h^. 

hi^nPh'^72''Rf,^p„ = -hi^h'^nPVpKp^ + K^pRPg + a^a^ 
-KhlV^cip. (A13) 

Equations (|XTT|) and (|XT3)) express the 4D curva- 

ture tensor in terms of projected quantities for a projec- 
tion operator based on an arbitrary spacelike vector 
When is hypersurface orthogonal, these expressions 
reduce to the usual Gauss-Codazzi equations. The case 
where is a KV is addressed in the next section. 

2. Field Equations expressed on the 3D quotient 
manifold in the case where 5'' is a Killing Vector 

We now specialize the results of the Gauss-Codazzi 
equations derived in the previous section to the case 
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where ^'^ is a KV obeying Eq. (pyTj) . When ^'^ is a KV 
the derivation and results presented here are equivalent 
to that found in [iOl. 

The Killing equation (12.11) implies that the tensor Kai3 
defined in Eq. (|A3[) is antisymmetric and thus the gener- 
alized Gauss equation (jA7|) holds. The Killing equations 
along with the identity $^A_^ = can be used to express 
the covector a/j defined below Eq. (jAlOp as 



(AM) 



This result, in conjunction with Eq. (jAlOp . then allows 
us to write the gradient of the KV as 



V^C. = -Ai/2i^^^--e[^V,]A. (A15) 



Substituting Eq. (jAlsp into the definition of the twist, 
Eq. (|2.2p . we obtain the expression 



(A16) 



which can be inverted using the identity e^i/po-e''^-^'^ = 
~6^[i/^P^CT] ^^'^ ^^'^^ ^^K^i, = 0, to yield Kpa in 
terms of the twist 



1 



2A3/2 



(A17) 



Finally substituting Eq. (|A17|) back into Eq. (|A15I) yields 
the identity 

V^e. = ^ W-C^w" - ^^[/.V.] A . (A18) 

Given the expressions for Kap, and V^^i/ in terms 
of the twist and norm of the KV, we can begin to evalu- 
ate the Gauss-Codazzi equations. An expression for Ricci 
tensor on the three manifold can be found by substitut- 
ing the double contraction of the KV with the 4D Rie- 
mann tensor, Eq. (|A13I) . into the contracted Gauss equa- 
tion (jA8[) and using the relations given in this section to 
arrive at 



1 



1 



1 

2"A 



4A2' 



This equation used extensively in Sec. Ill Al where it is 
referred to as Eq. (|2.7p . The first term on the right hand 
side of (|A13p vanishes because of the symmetry in the 
indices a, 5 on the left hand side of (|A13p and the anti- 
symmetry of K^a- 

Since is antisymmetric, K = Q. Substituting 

Eq. (|A17P and subsequently Eq. (jAlSp into the con- 
tracted Codazzi equation (|A12p we obtain 



(A19) 



which can be rewritten to yield 

which relates the twist of uj^ to the 4D Ricci curvature. 
This equations is reference as Eq. ()2.10p in Sec. Ill Al . 

The divergence of Wa is found by considering the to- 
tally antisymmetric part of the generalized Codazzi equa- 
tion, or equivalcntly contracting ^i.e""'''* with Eq. (jAlip . 
which becomes 



^pe^^'^^'n'^Ra^.^s = 2^^e^''^'{D5K^^ + a^K^s)- (A20) 

The first Bianchi identity, iJ^^po- -I- Rppua ~- Rg p^p = 
sets the term on the left hand side of Eq. (jA20p to zero. 
The right hand side of Eq. fJ20ll can be evaluated using 
the following expression for the derivative of the extrinsic 
curvature, 



D^K, 



"^^D^oj,. (A21) 



The resulting expression for the divergence of the twist 
vector quoted in Eq. (|2.9p is 



An equation governing the harmonic operator applied to 
A can be obtained by contracting the final projection of 
the Riemann tensor in Eq. (jAl3P with the three metric 
h"^ , and making use of the antisymmetry of and the 
expressions (|A17p for Kap and (|A14p for ap. The result 
is quoted in Eq. (j2.8p and given below, 

D^X ^^D^XD^X - \lo^lo^^ - 2R^,eC ■ 



2A 



A 



This completes the derivation of the reduced field equa- 
tions on S, used in Sec.|TT]to discuss axisymmetric space- 
times. 



Appendix B: Lorentz transforms of the 3D tetrad 

Here we discuss the effect of Lorentz transforms of the 
triad. As usual, these come in three types: boosts along 
the null vector P, and rotations about each of the two 
null vectors P and n^. First we discuss boosts. Let 



r = AP 



(Bl) 



Then li-j = Ali-j -f Ajli and hi-j = rii-j/A — Ajrii/A^. 
The nine rotation coefficients become 

7 7, 



e = e, 

e = Ae, 

a = Aa — A I 



A ^L, 
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(B2) 



where the directional derivatives are with respect to the 
original triad. The six curvature scalars transform as 



^5 = ^^05, 



00 ■ 



(B3) 



Next, let us consider rotations about the null vectors. 
The usual rotations about the null vectors in a null tetrad 
in the 4D are restricted to those that do not mix the axial 
KV with the three vectors that span the 3D hyper- 
surfaces which correspond to S. These transforms must 
leave the difference and sum of the complex NP spatial 
vectors to^ and m*'' invariant (since these correspond to 
the normalized KV and the spatial vector c'^, respec- 
tively). The usual rotations by complex parameters (see 
e.g. [261 ) are reduced to rotations by real parameters, a 
and b. 

For a rotation about we have 



r = V 



aP 



■n! + ac' + —r 



(B4) 



For a rotation about n* we have in complete analogy 



3^2 ^3 

bi = (j)i+ a[(j)Q + 04) -I- —03 + y 05 , 



00 = 00 + 2a03 -I- a05 . 



(B7) 



Finally, we consider rotations about n*. We first note 
that when interchanging V" and n*, the rotation coeffi- 
cients exchange identities as 



{a, /3, 7, 6, e, C, V, ^, 
-> {-S, ~C, -e, ~a, -(3, -r], -t, -6} . 



(B8) 



Thus, the effect of a rotation around n* by a factor b 
on the rotation coefficients can be derived from the ex- 
pressions given for a rotation around P by first applying 
the above relations to those transforms, and then tak- 
ing a b and swapping directional derivatives in the P 
and 71* directions, /^i — )■ /,2 and vice versa. The rotation 
coefficients thus transform as 



d = 5-bC, 

e^e-bS+-C-b^2, 
7~] = + b{S - l) -b'^C: 



C = C, 

^ = -/ + bL + b'^C, 

L = L + bC , 



e=n\ c'=c' + bn\ P^P + bc' + —n'. (B5) 

We are primarily interested in a fixed null direction P, so 
let us consider rotations about this vector. We have the 
following transforms for the rotation coefficients. 



a — ar 



7 = 7 — aa + —fi + aj 



2 ^ ' 



2 3 

J = 5 + a(r;-6) + y(c.-20)-y/3, 

2 3 
Qj CI 

1= L + a(7 - 77) + —{6 - 2a) + —(3 + 0,3 + aa^i , 



a 



+ a;2 + aa,3 + y",i 



(B6) 



In these expressions, the directional derivatives are with 
respect to the original triad vectors. The six curvature 
scalars transform as 



52 b^ 
a=a + 6(77 - 7) + -j{6 - 2i) - -C , 

^ =0 + 6(e - 77) + y (. - 25) + yC - & 3 - 

/3 + b{e - a) + y (7 + e - 277) + y (. - 

+ y C - bb^3 - y5,2 . 



66. 



(B9) 



We can write similar transformations for the cur- 
vature scalars 0j by noting that, under the ex- 
change of P and the curvature scalars transform as 

{05, 04, 03, 02, 01, 0o} -> {02, 04, 01, 05, 03, 0o} aild 

applying these transforms and a — >■ 6 to Eqs. (|B7p . 



Appendix C: Minkowski spacetime with as a 
coordinate 



To gain a better intuition into the choice of i/^ as a 
coordinate, let us consider Minkowski space in cylindrical 
coordinates (i, 0). In these coordinates p = and 
so dp = e^dip. Inserting this gives the line element 



ds^ = -dt^ + dz^ + e^'^id^p'^ + d(j)^) 



(CI) 



H + a03 + y 05 : 
h + a05 , 



Next, consider the metric on the conformal space S. We 
have 



% =diag[-e2^,e2^,e4'/'] 



(C2) 
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We can see immediately that R = 2e'^'^ , which trivially 
obeys the axis condition in 5 as "0 — )■ — oo. Looking at 
how the metric functions enter the triad in Eqs. ()5.15|) 
and (j5.16p . we see that an appropriate choice for a null 
triad adapted to the timelikc gradient Ta is 



r = (e-^e-'^,0)/V2, 
= (0,0,6-2'/'). 



n'-(e-^~e-^0)/V2, 
(C3) 



Using this triad, we can compute the rotation coefficients 
and begin to get a sense of the way each coefficient should 
behave as we approach the axis. However, first let us note 
that the triad chosen above has some troubling features. 
Comparing these to the corresponding tetrad vectors in 
M, as given by Eq. (|4.10p . we see that 

p = (e~'^,e-''',0,0)/V2, = {e^\-e^,0,0)/V2, 

= (0,0,6-^^,0). (C4) 

Near the axis, we see that our chosen and h'^ vec- 
tors are poorly behaved; l'^ blows up on the axis, and 
n'^ vanishes. We must boost the triad vectors by a fac- 



tor of A 



in order for the corresponding physical 



tetrad to the be well behaved on the axis. Computing 
the rotation coefficients on S with the boosted triad legs 
I' = (l,l,0)/\/2 and ^ e-2'^(l, -1, 0)/v^ leads to 



7 



(C5) 



with all others vanishing. Note that under a null boost, 
7 and e do not change; we cannot prevent the patho- 
logical behavior of these coefficients on 5 as -0 — oo. 
Meanwhile, jy does transform, and using a boost A = e™"^ 
we find fj = Tj ~- me^'^'^ , which can be used in this case 
to make rj — with the choice m ~ —1. This infinite 
boost at the axis has no effect on the vanishing of the 
other coefficients, and returns us to the triad we origi- 
nally considered in Ec^. (|C3p . 



Appendix D: Rotation coefficients associated with a 
triad adapted to the tp coordinate 

The rotation coefficients associated with the triads in 
Eqs. (|5.f3p (with h^p = 0) and (|5.16p and the corre- 
sponding metric Eq. (|5.f 5p are now computed using equa- 
tions (|4.5p and (j4.6p . Recall that for the case under con- 
sideration r] = ^(/3— 7+C~£) ^iid 7 = — e. The remaining 
coefficients are 



a 



V^ghs 

2ghs 
ikRs + hsRt) 
2V2gh,R ' 

_VR{ghs)^ 
2V2ghs 



\f2ghs 
\fR{hsnt^^ - nths.^) 

2ghs 
{utR^s - hsR^t) 



2V2ghsR 



(Dl) 



From Eq. (|D1[) the dominant scaling of the rotation co- 
efficients near the axis can now be obtained by requiring 
that the physical metric expressed in terms of {t, s, A, 0) 
coordinates is regular as the axis is approached, A — ?> 0. 

To examine the behavior of the metric components as 
we near the axis, we will quote a result of Rinne and 
Stewart [t^I- Consider a local Lorentz frame in a neigh- 
borhood near a point on the axis, p g and let us use 
Cartesian coordinates (x, y) on the space orthogonal to 
W2, so that the KV can be represented as 



-ydx + xdy 



(D2) 



If we insist that scalar quantities have a regular expansion 
in (x, y) about the axis. 



/(x,y)- E / 

m,n— 



(m, n)^m^n 



(D3) 



and that their Lie derivative with respect to vanishes, 
then it can be shown that the expansion must in fact be 
of the form 



J2 /^"^^" 

ri=0 



(D4) 



noting that X = + y'^ near the axis if our coordinates 
are appropriately normalized. By applying the same Lie 
derivative argument to the metric, it can be shown [72j 
that the {t, s) block of the metric admits expansions in A 
as if the metric functions are scalar quantities, 



9 s 



gts 



gtt 



n=0 



n=0 



(D5) 



It is always possible to choose s and t coordinates on the 
axis to be orthogonal to each other. This choice sets the 
first term in the off-diagonal metric function g^^"^ to zero. 

The series expansions about the axis in Eq. ()D5|) also 
set the series expansion for the functions entering into 
the metric hij in Eq. (j5.I5p . Explicitly we have that the 
functions hg, nt and It admit the following expansions 
near the symmetry axis: 



nt 



6^Enr)(t,.)A", Z,=6^Ezr^(t,.)A", 



n=0 



ri=0 



K = e^' E /ii")(t, s)A" g^e^'J] g^^Ht, s)A", (D6) 



where the coefficients p^"^ 



An) 



)/V2. If the met- 



ric is chosen to be diagonal on the axis, we further have 



that 1} 



(0) _ JO) 



In Eq. (|5.I9p we integrate the equations 



describing the rotation coefficient e by means of one of 
the Bianchi identities to yield {hsgY = g{t,s)/R. This 
result remains valid in the neighborhood of the axis and 
allows us to find an expression for the function g{t, s) in 
terms of the expansion coefficients h^^"^ and g^^'^ , namely 



2{h^)g 



(0))2 



(D7) 



33 



This result, together with the expansions given in 
Eq. (|D6p and the series expansion of R inipHed by the 
axis condition, 



(D8) 



allows us to determine that on the axis the rotation co- 
efficients scale as 



a 



y/g{t,s)y/X 



o(Vx), 



^ , ^ + O ( X^'A , 



O(A^), 



V5(M) 

C ^ ^ F== + O (Ai) , 



While (3 and C have O (A°) terms, we know that null 
rays which remain on the axis must be geodesic, and 
thus these terms must vanish. We then have 



(1) 



(1),(0) 



,(1) 



AO) 



AO) 



Q 



(0) 



(DIO) 



Recall that when working with a triad adapted to ip as 
a coordinate, Eq. (|5.7p states that the Weyl scalar ^'o = 
y/R/2(3, and as a result (jPlOp implies that the geodesies 
along the axis are principle null. It also implies that if 
the metric is chosen to be diagonal on the axis so that 



(0) 



J — n^t\ then this property persists to order O(A^), 

7(1) (1) 
smce /( = ni . 

We now substitute the expansions into the field equa- 
tions (|5.20p and (|5.2ip , and begin to solve them order by 
order in A. We start by looking at the subset of equa- 
tions that have directional derivatives only in the n° 
directions, namely Eqs. (|5.20p . and choose to set the met- 
ric diagonal on the axis. The dominant terms that arise 
from the sum and difference of the first two equations in 
Eqs. (|5?20l) give 



4/1 



(1)^ 



0, 



4/1, 



(1) 



,(0) 



= 0, 



respectively. These imply that 
.(1) 



WO)- 4^ 



(Dll) 



with fci a constant. Together Eqs. ((D9)) . (|DT0|) and ((DTT|) 
give 



e = -A-i - 4fci + 0{V\) 



(D12) 



The dominant A^^ term in the fourth equation of (|5.20p 
and the first equation in (|5.2ip governing 6,3 can only 
vanish if both ki — Q and the equation 



2h 



(0) 



-h2i?(i) = (D13) 



holds. With this, the equations in (|5.2ip that govern 
1^,3) 0^,3 J ^.3 1 s-nd t.3 are satisfied to 0{\fX) and the third 
equation in ((5?20)) to 0(A). Setting the 0(A°) term to 
zero in the equation governing 6 3 and the 0{\^) term to 
zero in those for C.3, /3^3 in (|5.2ip . and furthermore setting 
the 0(A°) term to zero in the fourth equation of (|5.20p 

fixes the R^'^\ np', if'^ and h')^^ coefficients in terms of 

h''s'\ R^^\ and their derivatives. These expressions 
arc lengthy, and we will only give the off diagonal term 
here. 



,(2) 



,(2) 



R 



(1) 

,ts 



16h, 



(0) 



(D14) 



Equation (|D14p , indicates that for a time-dependent met- 
ric the off-diagonal term It — rit is 0{X^^^). It is also clear 
that when the spacetime is dynamic, the diagonalization 
of the (t, s) block cannot be maintained off the axis. This 
off-diagonal term can be interpreted as a shift governing 
the motion of the coordinates along constant -0 slices as 
time progresses. Continuing on, it appears that the ex- 
pansion coefficients of the metric functions are fixed at 
each higher order by the lowest order terms , h^s\ 
and The same behavior occurs far from gravitating 

sources in asymptotically flat spacetimes, where the ex- 
pansion is in orders of inverse affinc distance. This is the 
asymptotic expansion of the Bondi formalism, which we 
discuss in Sec. IV Cl and Appendix IF] 



Appendix E: Geodesic null coordinates afflnely 
parametrized with respect to the 4D manifold 

In this appendix we write down the field equations 
adapted to a geodesic null coordinate system that is 
affinely parametrized with respect to the physical mani- 
fold A^. The results obtained in this section can be di- 
rectly derived from Sec. lVBl bv applying a boost A = e^^ 



using Eqs. (jB2p . and then selecting a new parameter 
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along the geodesic. However, because this choice of coor- 
dinates is used in our discussion of the Bondi expansion, 
we will give the equations in full here. Let us once again 
choose a null coordinate u such that h^^ujuj = 0, but 
this time we choose the null vector k ~ —e^^'u^i. The 
symmetry of u,a\b = ~(e^^'''/a)|(,, when expressed on the 
triad basis, leads to the following conditions on the rota- 
tion coefficients. 



(El) 



Making use of Eqs. (|4.12p . we immediately see that this 
choice for yields k = e = 0, so we see that is geodesic 
and affinely parametrized in the physical space. Let r 
denote the affine parameter in Ai, and choose the triad 
to be 

r = (0,1,0), n' = (e-2V'^_y;^^o), = (0, u;2, e""^')- 

(E2) 



Then the three metric on S becomes 

-e2* 



(E3) 



Applying the commutation relations to fmd r re- 

spectively yield 



7 
a 



-e. 



-^3,1, 



/3 = 0, 



ri = -2V',3 - e, 



(E4) 



and 



5 = -2wi?/>_i - wi^i, e = - {w29 - W2^i - 2^/1^3) , 

C = W2W3,2 + W2,2 + 2wiV',3 + ^1,3 . (E5) 

Substituting the expressions for the coefficients 7. a, /3, 
and 77 found from these relations into the field equa- 
tions (|3.15p -" P.23p and reorganizing gives 

1 =26iV',i + 2V'a + 6'^ (E6) 
e,i =27/.,iV^,3, (E7) 
(5.1 - 2V^,i (^-,2 + 2(5) - 2V'.i2 - (E8) 
561 - 6(4 - - (E9) 
267/.,! - (56*, (ElO) 

(Ell) 



f ,2 — e,3 
M + 02 
Ci + e,2 



-4CV',l +27^,2 V',3, 



C3 - ^2 = - 4C7^^3 - (5/, + ^2 , (E12) 
C,i + (5,3 + 2e,2 = Qe- 2it - 4C7/;,i - 2(57/>,3 

- 47^,2V',3 - 27/>,23 • (E13) 



As expected, the hierarchy present in Eqs. (|5.35p - (|5.42p 
(where 1°" is affine in 5) persists, which allows us to for- 
mally integrate the field equations. 

We conclude this appendix by detailing how the met- 
ric functions in Eq. ()E3I) transform with the remaining 



coordinate freedom [103j . These transformations include 
shifting the the origin of the affine parameter along each 
geodesic separately, t' = r + f{u,x), which transforms 
the metric function of (|5.27p according to 

w[ =wi- e^^'^f^u , W2 = "^2 + /xe~'"^ , W3 = W3 ■ 

(E14) 

Relabeling the individual geodesies withing a spatial 
slice, x' = diu, x) transforms the metric functions of (jE3p 
to 

/ I u)ti— 2)/i 9.U 1 2w^—M ( 9.u\ 

Wi — wi + e ^7f2— e * — 

9,x 2 Vs.x/ 

w'2=W2- e'"-^-2'/'^, e'"^ = (E15) 
9.x 9,x 

And finally, by relabeling the null hypersurfaccs, setting 
u' ~ h(u), t' ~ T/h,u ,the metric components transform 
as 



(E16) 



Appendix F: Derivation of the Asymptotic 
expansion of an affinely parametrized metric in null 
coordinates 

In this appendix we systematically solve the field equa- 
tions given in Appendix [E] in the asymptotic regime far 
from an isolated, gravitating system. Our method of so- 
lution further illustrates the integration of the hierarchy 
of the field equations that results when they are expressed 
on a null slicing. We focus o nly on spacetimes that ad- 
mit the peehng property [S^-ISa]. In this case the Weyl 
scalars have a power series expansion at future null infin- 



ity of the form = t'-^ Y.n=o where ^"f"' are 

constant along an out-going null geodesic, i.e. '^^/^\u, x), 
and T is the affine parameter along the geodesic. Using 
Eq. (|5.46p as a starting point, we derive the power se- 
ries expansion of the metric functions and the rotation 
coefficients in terms of a series in 1/r. This information 
makes apparent the boundary conditions that have to 
be imposed when solving the complete set of equations, 
and provides explicit information about the fall-off of all 
rotation coefficients at null infinity . 

Consider first the properties of the null tetrad vector 
Z^. We will take to be the tangent to out-going, null 
geodesies far from the isolated system, so that (3 — 0, and 
also affinely parametrized in A4, which sets a = —2'ip,i. 
Directional derivatives in the P direction on S can thus 
be expressed as f^i = f^r- The simplified field equations 
and form of the metric for this choice of parametrization 
is given in Appendix [E] Specifying the series expansion 
for on a null hypersurface of constant u gives almost 
all of the data required to continue the spacetime off the 
hypersurface. During the calculation that follows, we 
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quantify how this information is transmitted to the met- 
ric functions. 

From Eq. ()4.15|) and (j5.46|) we have that the leading 
order terms of the function ip are related to the coefficient 
via + {ip,if = r-s^-^,"^ + 0(r-6). Solving this 
equation term by term, we find that 

(Fl) 

Here, 

^(0) ^ ^(0) 

(m,x) is a function whose properties 
have yet to be defined. We will see that cr^°'' corresponds 
to the dominant term in the series expansion of the shear 
of i^. Its labeling corresponds to the choice made in jl23l | , 
whose derivation we initially follow closely when working 
out the expansion properties of the optical scalars. In 
the scries expansions that follow we keep terms of suffi- 
ciently high order to indicate where dominant terms of 
the expansions of the Weyl scalars enter into the metric 
functions. 

Integrating Eq. (|Fip with respect to r adds an addi- 
tional integration constant 'i/''-'^-' (u, x), and allows us to 
express ?/; as 



(0) 



ij} =i'^^^ +\nT 

T 

+ 0(r-4). 



^(0)2 (^(0)3 _ ^(0)^2) 



2r2 



3r3 



(F2) 



A series expansion for 6 can be found by substituting a 
power series ansatz for 6 into Eq. (|E6|) and using Eq. (jFl|) 
to define the series expansion for -0.1. It can be shown 
that the leading order behavior of the solution admits 
only two possibilities, namely = —t~^ + 0{t~^) or 
9 = — 2r~^ -|-0(r~^). The former corresponds to a cylin- 
drical type spacetime that is not astrophysically relevant. 
We will only consider the latter case and further discuss 
the justification for this choice after Eq. (jF4p . Using our 
coordinate freedom to relabel the origin of the affine pa- 
rameter, t' = T + f{u, x), it is always pos sible to set the 
next coefficient in the expansion to zero jl23| . Examin- 
ing the remaining coefficients in Eq. (jE6|) term by term 
leads to the scries expansion 

, (.^..y^^ 

With this, we can integrate Eq. (jE4|) . w^^i = —9, to arrive 
at an expression for the metric function ws, 

W3 =w'i'> + 2 ln(T) 2 



6r4 



+ 0(t 



-5\ 



(F4) 



which adds the integration constant w^\u,x) to our 
list of undetermined expansion coefficients. Note 
that the leading order term of the metric coefficient 
9xx ~ e^™^"^''' in A4 is proportional to r^, which is typ- 
ical for a surface of constant r that is asymptotically 



spherical. Had we made the selection 9 = —t~^ +0{t~'^) 
above, the leading order term would have been indepen- 
dent of T. 

The next set of variables to be considered in the inte- 
gration hierarchy are e, ^i, and the metric function W2- 
For our chosen triad, determining e also fixes two other 
rotation coefficients; from Eqs. (|EI[) and (jE4|) we have 
that r] — —2-0^3 — e and 7 = — e. The peeling property of 
the Weyl scalars (j5.46p in conjunction with the expres- 
sion for ^1 given in Eq. (|4.I5p and the field equation (|E7p 
yields 



2e0,; 



4e-'^ 



V2 



[^fK-^ + Oir-')), (F5) 



which can be systematically solved to yield a series solu- 
tion for e, 

g(0) g(2) g(3) 



Next, using Eq. (|E7[) . we find a series expansion for ip^^ 
of the form 



c(o) ^(0)£(0) _ 2£(i) 2(7(0)e(i) - 3e(2) 



2t 



2t2 



6^(0),(2)_ ^(0)^(0) _ 8^(3) 

47^ + ^ 



2r3 
_-5^ 



(F6) 



The higher order coefficients in the expansion for e are 
fixed in terms of existing quantities as follows. 



(3e(i 



e(o)^(o) 



6 



3gV'(o) ' 
(F7) 



and so far the coefficients e^*'^ and e^^'' are unconstrained. 
It turns out that the leading coefficient e^^^ can be set to 
zero using a gauge transform. 

To see this, note that we can obtain an expansion for 
the metric function W2 using Eq. (|E5|) . The resulting 
expansion is 



,(0) 



(0)^(0) 



W2 



~2 



w. 



(2) 



6r3 



+ 0(r-4). 
(F8) 



It is then possible to make use of the gauge transfor- 
mation x' = g{u, x) to relabel the geodesies such that 
6^°' = 0; using Eq. (jE15[) . The resulting expansions for 
e. -0^3 and W2 reduce to 



.(1) 



I 



2aW 3a(0)^ 



(0) 



(F9) 
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W2 = 



(0) 



„(2) 



4 



3^V^ 



+ 0(r-4). 



(FIO) 
(Fll) 



The next set of variables obtained via the systematic 
integration of the field equations includes (5, the Weyl 
scalar ^'2, the metric function wi, and the derivative 
■0^2 ■ The coefficient 5 can be obtained by integrating 
the field equation (|E8|) . where the commutation relations, 
Eq. p.29|) . the definition (|4.15p of #2 and the harmonic 
equation p.30p for ij), have been used to replace direc- 
tional derivatives in the direction with known series 
expansions. The resulting equation 



1 + 2(SV,i = (e + Vas)^ - 2e-2'A^2 



(F12) 



implies that 5 must admit the series expansion 

,5(1) 2a("'<5(i) 3aW2 5(i)+e-2V''°'^(0) 



(F13) 



The expansion for 5 can now be used to obtain the 
expansions for the metric function wi using Eq. (|E5p . 
i5 = —2'wi'ip^i — wi.i. The first two terms in the resulting 
expression are 



,(2) 



Wi 



0{t- 



(F14) 



If the metric g^^ is to be asymptotically fiat, the metric 
function e^'^wi must be finite as r 00, and thus the 
integration constant (j'^' ~ 0. The resulting expansion 
for the metric function wi becomes 



Wi 



(0) 



0(r-4). (F15) 



A series for the directional derivative V',2 can be obtained 
from the field equation (|E8[) using the commutation re- 
lation to switch the order of differentiation on tp. The 
resulting expression becomes 



6V',iV',2 + 2i/),2i = —2Sil,\i — 5 1 



(F16) 



which implies that V-'.2 has the scries expansion 



-2^(«)^(0) 



„(2) 



(2a(°V 



(F18) 



Examining the directional derivative ip^2 and the expan- 
sion pTTl) . we obtain ^(i) = -e-2'>"'V^°\« = which 
implies that i/'^^) = '(/''•^•'(x) is independent of u. We also 
have that 



(0) 



-w 



(2) _ 2^(«)^(0) 



(F19) 



Finally, we examine the remaining field equations to 
obtain further restrictions on the expansion coefficients. 
Substituting the expansions and directional derivatives 
obtained thus far into Eq. (|E9[) yields the condition 



(0) 



27/;(o)_ 



0, at order r ^, which implies that 



w- 



(0) 



(0)/ 



(F20) 



Note however that the coordinate transformation 
x' ~ 9ix) transforms the metric function as e^^ = 
e"'''/5.x and allows us to set /(°' (x) to any arbitrary func- 
tion of our choosing. We can understand the meaning of 
a choice of /C^) by insisting that as r — ?> 00, the (x, 0) 
block of the physical metric on Ai has the geometry of a 
sphere. In other words, qab — >■ ^ab, where ^Iab is a 
metric on the unit 2-sphere, and {A,B} € (%,?!>) • With 
this requirement on the angular geometry of out-going 
null surfaces, a particular choice of /(°) allows us to fix 
both ■f/''-^) and to identify the particular angular coordi- 
nate X corresponding to this choice of /'■'^•'. For example, 
setting /'■^•'(x) = — 2?/'(°' corresponds to the choice where 
the metric on the unit 2-sphere ^Iab has unit determi- 
nant. In this case, we find that e^''^^"^ ~ 1 — x^i where 
X = cos 9, 9 is the usual angular coordinate, and the axis 
is located at x = il- For comparison, note that setting 
f^^Hx) = corresponds to the Fubini study metric repre- 
sentation of the sphere used in |l23l | . Henceforth we will 
set /(o)(x) = -2V'(°\ which means that we have selected 
X = cos 9 and that w^^ = 0. 

Evaluating Eq. (jE9p at order yields the additional 
condition 



.(2) 



_^|0)2_1^(0)^ 



XX 



2(1 



(F21) 



Next, at order r ^ we obtain an equation that evolves 
,(0) 



0(t-^). (F17) ^-J,"^ from one hypersurface to the next hypersurface 



Since the series expansions for ip and the three metric 
functions wi, W2, are given, the directional deriva- 
tives of any function expressed as a series in 1 /r can also 
be expanded as a series. We begin by examining the di- 
rectional derivatives of ip to determine what restrictions 
the resulting expressions place on the existing expansion 
coefficients. By considering the directional derivative ■0.3 
and the expansion (|F10p we obtain the result 



.(1) 



(0) 



V2 



(F22) 



Using equations (|E4[) and (jE5|) . the expansions for l and 
C can be shown to have the form 



0\ 



,,(0) 



,,(1) 



(F23) 
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where the expansion coefficients are related to those func- 
tions already defined by 



,(0) 



1 



(0) 
2 



(F24) 



Using Eq. (j4.15p the dominant terms in the remaining 
Weyl scalars can be shown to be 



(0) 



v]/ 

^3 ^2 



+ 0(r-2), 



where 



(0) 



[(i-xV"\J 



2^1 -X^ 



(F25) 



(F26) 



The evolution equations that propagate the coefficients 
and from one null hypersurface to the next can 
be obtained by examining (jElip and (jE12p respectively 
at 0(t~^), yielding the expressions 



(0) 
2 ,X 



2^1 -X^ 



(0) 



(F27) 



(F28) 



The field equations (|Ef Op and (|Ef 3p yield no additional 
constraints and vanish to 0{t^^). Also note that the 
harmonic equation for ip has been satisfied. 

The results obtained thus far are now briefly summa- 
rized. The 4D line element can be expressed as 



ds^ = ~2e^'^'wi du^ - 2du dr + e"'^W2 du dx 
+ e^-^-^-^v. ^^2 _^ g2^ ^^2 ^ 



where, according to Eqs. (|F4|, ((FTT|) . ((FT5|) . ((FT8|) . 

and (|F2f [) . the metric functions admit the following 
asymptotic expansion as the affine parameter t — > oo: 



„2V 



(0)' 



(f - x') ( - 2raW + + ^ ) + ^(r-^), 



Wi = 



W2 



f 

2^2 



6t2 



(0) 



(1 - X^) 

[(i-x^V^°)] 

(f-X2)T2 



+ 0(r-^), 



2V2*i 



(0) 



3^1 



X" 



+ 0(r-4). 



(F30) 

To fully specify the solution we must make a choice for 
the functions ^!f'\u^ x): * G {07 1:2} on a hypersurface 
and specify the function cr'^°) (w, x) for all w and x- A 
natural hypersurface to choose is u — >■ cxj and to specify 
the functions to correspond to the Schwarzschild solution. 
In this case 



{l'o,^'i,*2,<7(°',<7(°)} ^ {0,0, -M, Ce 



-2^(' 



,0}. 
(F31) 



The constant M is the mass of the final black hole and C 
is an arbitrary constant. The fact that the shear of the 
null bundle must be regular on the axis sets C = 0. For 
u < oo, ^l*'^ can now be determined provided a^^^{u, x) 
is given. 

The results given this section are further discussed in 
Sec. ESI 

The expansion of the 4D spin coefficients can be ob- 
tained from the results in this appendix using Eqs. (j4.12p . 
The series expansion of the shear a is found to be 

(0) (0) 3 _ ^(0) 

^- — + ^- + 0{T-% (F32) 



confirming that is indeed the dominant term in the 
expansion of the shear. It should be noted that our ex- 
pansion for the spin coefficient f is 



-/3 = -^(VA3 + e) = 



0{t- 



(F33) 



and that the prefactor of 1/3 in front of the r ^ term 
differs from the result obtained in Il23t. 
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